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ABSTRACT 

We present the results from a survey for 12 CO emission in 40 luminous sub-millimetre 
galaxies (SMGs), with 850-/im fluxes of Ssso^m =4—20 mJy, conducted with the 
Plateau de Bure Interferometer. We detect 12 CO emission in 32 SMGs at z ~ 1.2 - 
4.1, including 16 SMGs not previously published. Using multiple 12 CO line (J up = 2- 
7) observations, we derive a median spectral line energy distribution for luminous 
SMGs and use this to estimate a mean gas mass of (5.3 ± 1.0) x 10 10 M . We report 
the discovery of a fundamental relationship between 12 CO FWHM and 12 CO line 
luminosity in high-redshift starbursts, which we interpret as a natural consequence 
of the baryon-dominated dynamics within the regions probed by our observations. 
We use far-infrared luminosities to assess the star-formation efficiency in our SMGs, 
finding a steepening of the L^q-Lfir relation as a function of increasing 12 CO J up 
transition. We derive dynamical masses and molecular gas masses, and use these to 
determine the redshift evolution of the gas content of SMGs, finding that they do not 
appear to be significantly more gas rich than less vigorously star-forming galaxies at 
high redshifts. Finally, we collate X-ray observations, and study the interdependence 
of gas and dynamical properties of SMGs with their AGN activity and supermassive 
black hole masses (Mbh)j finding that SMGs lie significantly below the local Mbh-ct 
relation. We conclude that SMGs represent a class of massive, gas-rich ultraluminous 
galaxies with somewhat heterogeneous properties, ranging from starbursting disc- like 
systems with L~ 10 12 Lq, to the most highly star-forming mergers in the Universe. 

Key words: cosmology: observations - galaxies: evolution - galaxies: formation - 
galaxies: ISM 



1 INTRODUCTION 

The discovery of a population of sub-millimetre (sub-mm) 
bright, highly star-forming galaxies at high redshift has pro- 
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vided a critical challenge for hierarchical galaxy evolution 
models (Blain et aL]|1999| |Baugh et al]|2005[). These lu- 



minous "sub-mm galaxies" (SMGs) are young, highly dust 
obscured galaxies, at a median < z >~ 2.5 (e.g. |Chapman 
20031 I Chapman et~aLl|2005l hereafter C05; |Wardl~ 



et al. 



et al.||2011 1, with extreme far-infrared (far-IR) luminosities 
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(Lfir = 10 12 -10 13 Lq) i mplying prodigious sta r formation 



rates of £ 10 3 M yr _1 ( Chapman et al.||2010 Magnelli et 
al. 2012). As they are selected at 850/xm (corresponding to 
restframe A ~ 200-400/im) these sources tend to have signif- 
icant masses of cold dust (C05, see also Magnelli et al. 2012), 
as well as correspon dingly substantial reservoirs of molecu- 
10 " 



lar gas (~ 10 iviq; 
high stellar masses ( 



Greve et al.|2005| |. Together with their 

11 mV 



Hainline et al. 



20111 



10" Mq; e.g. 

this suggests this population reside in some of most massive 
dark matter haloes in the high-redshift Universe. 

SMGs have been studied extensively since their discov- 
ery in the 850-/xm atmospheric window using the Sub-mm 
Common User Bo lometer Array (SCUBA) on the JCMT 
( Smail et al.|1997 \ , and it is clear that they represent a pop- 
ulation of cosmological importance. Recent advances in IR- 
detector technology have opened up the possibility of much 
wider field surveys at these wavelengths, including Herschel 
SPIRE at 500^m; (Chapman et al. 2010; Magnelli et al. 
2010, 2012), in the millimetre- waveband with the South Pole 
Telescope (Vieira et al. 20101 and surveys with the new 



SCUBA-2 camera (Holland et al. 2006) on JCMT. These 
surveys are providing large numbers of dusty galaxies with 
properties similar to the original SCUBA population. In this 
work, we use the term "SMC" to refer to 850-pim-selected 
galaxies which comprise the bulk of our sample. 

SMGs lie above the "main se quence" of star formation 
(denned as SFR/M* ~ M?(l+z^ 



Schiminovich et al. 2007 



|Noeske et al.|2007||Daddi et al.|2007||Rodighiero et al.|20l'o] 

SMGs typically lie ~ 1 dex above this "main sequence" ; sig- 
nificantly more than the scatter around the sequence, which 
is typically ~ ±0.3 dex. 

As a result of their location in the SFR-M* plane (along 
with corroborating kinematic and morphological evidence; 
Tacconi et al. 2006, 2008; Engel et al. 2010), many authors 
have argued that SMGs can be understood as "scaled-up" 
analogues of ultra-luminous IR galaxies (ULIRGs) in the 
local Universe, which typically lie an order of magnitude 
(or more) above the SFR-M* "main sequence" . In this pic- 
ture, luminous SMGs are best understood as major-merger 
events, with the extreme star formation being driven by the 
merger. 

There has also been claims, however, for a secular origin 
for some SMGs. Theoretical simulations that model SMGs 
as massive, star-forming discs reproduce some of the physical 
properties of the SMG population, while avoiding some of 
the possible difficulties that arise from ascribing a merger- 
origin to all SMGs (such as the observed number counts). 
While there has been some corroborating evidence for this 
picture (e.g. Bothwell et al. 2010 presented evidence for a 
luminous SMG with disc-like properties) , simulated disc-like 
SMGs generally fail to reproduce the extreme star formation 
rates characteristic of the bright end of the SMG population. 

In addition to these (apparently opposing) pictures, 
"main-sequence" galaxies can also fall within the stan- 
dard SMG selection criterion (a simple flux cut at 850/im, 
S850(jm ^ 4mJy). High-resolution imaging has resolved a 
few sub-mm sources into two (or more) less luminous star 
forming galaxies, which happen to fall within the same sub- 
mm detection beam, pushing the combined system above 
the SMG survey detection limits (see |Wang et al.||2011 |. In 
addition, Chapman et al. (2005) confirmed the presence of 
a population of lower luminosity galaxies at lower redshifts, 



which are selected at sub-mm wavelengths due to unusually 
cold dust temperatures. An overarching theme of the results 
over the last few years has therefore been an appreciation 
of the diversity of the SMG population. The transition be- 
tween the brightest starburst sources and the more secular 
systems is a undoubtably gradual one, and it is clear that a 
mix of models is required if we are to fully understand the 
entirety of the SMG population. 

Observations of the cold molecular phase of the ISM 
have the power to resolve many of these issues, as well as 
providing unique insight into the physics and behaviour of 
these systems. Molecular line observations provide the most 
direct insight into galaxy dynamics (e.g. Walter et al.|2011 \ 
with turbulent merging systems having very different kine- 



matic profiles from ordered rotating discs (Neri et al.|2003 



Tacconi et al.|2006| 2008, 2010). In addition, molecular gas 



linewidths provide excellent estimates of the dynamical mass 
of galaxies, free from the uncertainties (extinction effects 
and non-gravitational motions, such as winds) that can po- 
tentially affect estimates derived from optical and near-IR 
observations. 

From the intensity of the observed 12 CO emission lines, 
it is possible to calculate the mass of the molecular gas 
reservoir. As molecular hydrogen lacks a permanent elec- 
tric dipole, it is all but invisible for even nearby galaxies. 
Instead, 12 CO is used as a tracer molecule, as it is produced 
in similar conditions to molecular hydrogen and is collision- 
ally thermalised by H2 even at low densities (though it is 
typically optically thick). Once the relative ratio of 12 CO/H2 
molecules is known (the conversion factor is usually param- 
eterised as a), it is possible to calculate the mass of the 
underlying H2 reservoir. When combined with far-IR lumi- 
nosities (calculated directly, or inferred from radio contin- 
uum flux observations via the far-IR/radio correlation), gas 
masses provide estimates of the star formation efficiencies in 
SMG s (e.g. Greve et al. 2005; |Genzel et al.|2010||Daddi et"aT 



2010 1. Whether such extreme starburst galaxies have simi- 



lar star formation efficiencies to more "normal" star forming 
galaxies at similar epochs, or whether modification is needed 
in order to explain their prodigious star-formation rates, has 
important implications for galaxy evolution models. 

Historically, a number of instrumental limitations 
caused molecular gas in high-redshift galaxies to be poorly 
studied. Observing molecular gas is very time consuming 
with current instrumentation, with on-source integration 
times of several hours being typical for even luminous 12 CO 
sources. In addition, several steps are required in order to ob- 
serve a 12 CO line once an SMG has been detected at sub-mm 
wavelengths. After the galaxy has been detected with the 
typically low-resolution sub-mm or far-IR beam, an accu- 
rate position is required using radio continuum emission (Ivi- 
son et al. 2001; Biggs ct al. 20111. Only then can UV spec- 
troscopy obtain an accurate redshift (C05), and only at that 
point could the relatively narrow-bandwidth receivers ob- 
serve 12 CO transitions. Of course, with the higher-resolution 
imaging in the far-IR/sub-mm made available by the advent 
of ALMA, the positioning of SMG-bright sources is becom- 
ing easier, and with increased bandwidths available at mm 
wavelengths the required redshift accuracy is now becoming 
less of an issue). 

One significant weakness of the studies of the gas con- 
tent of high-redshift starbursts published to date is their 
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poor redshift coverage, and subsequent limited ability to 
address evolutionary effects in the molecular gas properties 
of SMGs. The redshift evolution of molecular gas fractions 
and dynamical masses encodes important information relat- 
ing to the assembly of massive galaxies, providing valuable 
constraints on models of structure formation. 

In view of the importance of this population, and of 
the paucity of existing datasets, we undertook a large sur- 
vey for 12 CO emission in a large sample of SMGs using the 
IRAM Plateau de Bure Interferometer (PdBI), aiming to 
draw significant conclusions about the SMG population as 
a whole across a wide range of redshifts. This makes it pos- 
sible, for the first time, to address evolutionary effects in the 
gas properties of luminous SMGs. 

Our IRAM PdBI SMG survey was started in 2002. 
Prior to the start of the survey, only two SMGs had been 
detected in 12 CO; SMM J02399— 0136 (z = 2.81: Ivison 



et al. 1998; Frayer et al. 1998 


Genzel et al. 2003 1 and 


SMM J14011+0252 (z = 2.56: 


Frayer et al. 1999 Ivison 



were published by Neri et al. ( 2003 \ and Greve et al. ( 2005 1, 
who added 12 CO detections of three and five new SMGs 
respectively. In addition, spatially resolved (sub-arcsecond) 
imaging for a total of 12 SMGs from the survey has been re- 



por ted by Tacconi et al . (2006, 2008), |Bothwell et al.| ( |2010[ ), 
and |Engel et~aT] p010| . The IRAM PdBI SMGs survey has 
now reached its conclusion with a total of 40 SMGs at z ~ 1- 
4 observed in a range of 12 CO transitions. 

In ij2]we give details of the survey, the observing strat- 
egy used, and the reduction and analysis of the data. S|3] 
presents a median spectral line energy distribution for the 
sample, from which luminosity ratios are estimated. In Q 
we describe the physical properties of the SMG sample 
as revealed by the 12 CO observations, including the kine- 
matic properties, dynamical masses, star formation efficien- 
cies and molecular gas properties. In §4.6| we also compare 
baryonic and dynamical mass measurements, and discuss 
the implications for deriving physically-motivated parame- 
ters for SMGs. SjSJthen discusses the effects of supermassive 
black hole (SMBH) activity on the SMG population, and !|6] 
presents our conclusions. 

Throughout this work we adopt a cosmological model 
with (Q m , fl A , H p) = (0.27, 0.73, 71 kms" 1 Mpc -1 ), and 
a 



Chabrier (2003) IMF 



2 OBSERVATIONS, REDUCTION AND 
SAMPLE PROPERTIES 

2.1 Sample Selection 

Our full sample of SMGs observed with the PdBI low- 
resolution programme consists of 40 galaxieij]] at z ~ 1— 
4. The initial target selection was drawn from the optical- 
spectroscopic survey of SMGs by Chapman et al. (2005), 
as well as the SCUBA Cluster Lens Survey (Smail et 



1 Four SMGs in our sample are members of 

close pairs (SMM J123711.98/SMM J123712.12, and 

SMM J123711.86/SMM J123708.80), and in each case a sin- 
gle pointing was sufficient to observe both SMGs. 



al. 2002) and the SHADES survey of the Subaru/XMM- 
Newton Deep Field (Coppin et al. 2006). Given the dura- 
tion of the survey, our observing strategy naturally evolved 
with time, taking advantage of the rapid advance of other 
multi-wavelength projects studying this population. There- 
fore as the survey progressed we also included a few of the 
first spectroscopically-identified millimetre-selected galaxies 
(from surveys with the MAMBO camera on the IRAM 30m; 
Bertoldi et al.|2000| [Greve et al. 12004] ), as well as SMGs 



with more precise rest-frame optical redshifts derived from 
near- or mid-IR spectroscopy. As a result the final tar- 
get sample includes SMGs selected to have precise spectro- 
scopic redshifts, derived using one or more of the following 
techniques: (i) Lyman-cv emission (C05); (ii) Ha emission 
( jSwinbank et al. 20041, and (iii) PAH emission from mid- 
IR spectroscopy with Spitzer ( Pope et al. 2008 ; Mcncndcz- 
Delmestre et al.|2009 i. 

For those SMGs with only optical (restframe UV) spec- 
troscopy at the time of observation, we used the median ve- 
locity offset between optical and near-IR redshifts for SMGs 
to estimate the likely velocity offsets of their systemic red- 
shift from the optically-derived value. The resulting redshift 
uncertainty is Az = 0.005 and this should guarantee that 
the 12 CO emission from an SMG at z ^ 1 would fall within 
the 580-MHz bandwidth available at 3 mm at the time. 

As we have radio fluxes for our sources, we can derive 
the far-IR luminosity for our SMGs from the 1.4-GHz con- 



tinuum, via the far-IR-radio correlation in Yun et al. ( 2001 1 

;.4 x io 14 ) Si.4 (i + 2) (a-1) , 



47rDf 



(1) 



where Dl is the luminosity distance in metres, Si. 4 is the 
radio flux density at 1.4 GHz (in Wm _2 Hz _1 ), and a is 
the synchrotron slope used to K-correct the 1.4 GHz obser- 
vations to the appropriate source redshift (a is taken here 
to be 0.8). This equation assumes 9fir = 2.34. Some re- 
cent results from Herschel have suggested that a lower value 
(?fir ~ 2.0) may be appropriate for luminous SMGs (Ivi- 
son et al. 2010; Magnelli et al. 2012), which would lower our 
derived far-IR luminosities by a factor of ~ 2. However, no 
clear conclusion has yet been reached on this issue, so fol- 
lowing Magnelli et al. (2012) we assume the (local) value of 
(/fir = 2.34, and note the possibility of our far-IR luminosi- 
ties being over-estimated. 

2.2 Data Acquisition, Reduction and Analysis 

The observations of SMGs described here were obtained in 
observing campaigns on the Plateau de Bure Interferom- 
eter between 2002 and 2011, in good to excellent weather. 
The observations were undertaken with the lowest-resolution 
"D" configuration of the interferometer, in order to max- 
imise the sensitivity, and used five of the six available an- 
tennae (giving a total of 10 baselines). The targets were 
observed in the 2- or 3-mm windows, depending upon the 
redshift of the source. The typical resulting resolution of our 
3-mm maps is ~ 4". Table 1 lists our observation log. 

For data reduction, the IRAM GILDAS software was 
used. Data were carefully monitored throughout the tracks, 
and any bad or high phase noise visibilities were flagged. 
One or more bright quasars were typically used for passband 
calibration, and flux calibration was determined using ob- 
servations of the main calibrators, 3C 454.3, 3C345, 3C273, 
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|03 02 27.763 +00 06 53.0o| ; 




30 20 10 -10 -20 -30 

RA (arcsec) 




-1000 -500 500 1000 

Velocity (km/s) 



30 20 10 -10 -20 -30 

RA (arcsec) 

Figure 1. Example maps and 12 CO spectra for two of the SMGs in our sample (the full sample is shown in Appendix A). Top: 
SMM J030227.73. Bottom: SMM J105141.31. The maps show 12 CO (blue contours) and 24-fim (red contours) overlaid on near-IR images 
made from the combined Spitzer IRAC images (3.6, 4.5, 5.8, and 8^tm). 12 CO contours start at 2-<r significance, and are spaced in steps 
of Ict (top) and 2<r (bottom). The cross shows the radio position, and the ellipse in the bottom right shows the synthcsised PdBI beam. 
The spectra are taken at the point of maximum 12 CO intensity. The best-fitting single Gaussian profile is overlaid on each. Vertical 
dashed lines indicate redshift estimates from the literature: red shows redshifts derived from nebular recombination lines (Ha, or 0[n]), 
black is from UV emission lines (C05). 



and MWC 349. Phase and amplitude variations within each 
track were calibrated out by interleaving reference observa- 
tions of nearby compact calibrators every twenty minutes 
(see Table 1). For the subsequent analysis, the routines CLIC 
and mapping were used, producing naturally-weighted dat- 
acubes which were then outputted for analysis with our own 
idl routines. 

In a few line was detected close to the edge of 

the band pass in the initial track. In these situations the 
frequency setting was then adjusted to centre the line in the 
bandpass, and the source was re-observed. On average dur- 
ing the survey, a source was typically observed for 2-3 tracks 
(a total on-source integration time of 10-18 hours). If no sig- 
nal had been detected after this, the source was classed as 
a non-detection. Of the non-detections, subsequent observa- 
tions (either in 12 CO, or from spectroscopy - either optical 
or near-/mid-IR) have shown that some of the originally 
targeted redshifts were incorrect and that the 12 CO lines 
likely fall outside of our band (Table [2|. We therefore con- 
clude we have only eight true non-detections. Table 4 lists 



the observed 12 CO properties of our SMGs. Figure [T] shows 
example maps and spectra for two of the SMGs in our sam- 
ple - the full sample of detections and non-detections are 
shown in Appendices A and B respectively. 
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Table 1. Observation log 



ID 








u obs 

(GHz) 


RMS 
(iii 


Banc] | 


RMSchannej^J 

(mjy) 


On-source time 
(hours) 


Calibrator 


Reference 


SMM J021738 


.71- 


-050339. 


.7 


151.9 





.15 


0.80 


4.1 


3C454.3 


This work 


SMM J021738 


.91- 


-050528. 


.4 


130.6 


0. 


.15 


0.80 


5.1 


3C454.3 


This work 


SMM J021725 


.16- 


-045934. 


.7 


139.5 


0. 


15 


0.80 


5.3 


3C454.3 


This work 


SMM J030227 


.73- 


h000653. 


.3 


95.73 


0. 


15 


0.80 


13.1 


0336-019 


This work 


SMM J044315 


.00- 


K021002. 





150 





.17 


0.7 






Neri et al. (2003) 


SMM J094304 


.08- 


h470016. 


2 


150 





.17 


0.7 






Neri et al. (2003) 


SMM J105141 


.31- 


H571952. 





104.12 


0. 


12 


0.68 


19.3 


0923+392 


Engel et al. (2010) 


SMM J105151 


.69- 


^572636. 





87.92 


0. 


05 


0.51 


8.7 


0954+658 


This work 


SMM J105227 


.58- 


K572512. 


.4 


99.94 





.09 


0.84 


7.1 


0954+658 


This work 


SMM J105230 


.73- 


h572209. 


.5 


96.02 





12 


0.65 


21.9 


1044+719 


This work 


SMM J105307 


.25- 


^572430. 


.9 


91.33 


0. 


05 


0.38 


16.4 


0954+658 


Bothwell et al. (2010) 


SMM J123549 


.44- 


h621536. 


.8 


107.95 


0. 


12 


0.65 


11.4 


1150+497 


Tacconi et al. (2006) 


SMM J123555 


.14- 


h620901. 


.7 


80.46 


0. 


.19 


1.0 


11.0 


1044+719 


This work 


SMMJ123600 


.10- 


h620253. 


5 


115.49 





.42 


2.2 


9.4 


1044+719 


This work 


SMM J123606 


.85- 


K621047. 


2 


98.65 


0. 


.14 


0.78 


18.0 


1150+497 


This work 


SMM J123618 


.87- 


i-621007. 


.5 


107.72 


0. 


06 


0.57 


13.0 


1418+546 


This work 


SMM J123618 


.33+621550. 


.5 


153.73 


0. 


.09 


0.64 


18.6 


1044+719 


Bothwell et al. (2010) 


SMM J123621 


.27- 


K621708. 


.4 


154.19 





15 


1.0 


7.3 


0954+658 


This work 


SMM J123629 


.13- 


K621045. 


.8 


114.52 





.34 


1.8 


6.1 


1150+497 


This work 


SMM J123632 


.61+620800. 


1 


115.48 


0. 


.54 


2.9 


20.1 


1044+719 


This work 


SMM J123634 


.51- 


h621241. 


.0 


103.71 


0. 


.13 


0.79 


15.1 


1044+719 


Engel et al. (2010) 


SMM J123707 


.21- 


h621408. 


1 


99.08 





16 


0.85 


9.6 


1150+497 


Tacconi et al. (2006) 


SMM J123708 


.80- 


K622202. 





159.86 





21 


1.0 


17.5 


1044+719 


This work 


SMM J123711 


.98+621325. 


.7 


115.41 


0. 


.33 


1.8 


12.3 


1150+497 


Bothwell et al. (2010) 


SMM J123711 


.86+622212. 


.6 


159.86 


0. 


21 


1.0 


17.5 


1044+719 


Casey et al. (2009) 


SMM J123712 


.05- 


K621212. 


3 


88.37 





08 


0.55 


17.8 


1044+719 


This work 


SMM J123712 


.12- 


K621322. 


2 


115.41 





.33 


1.8 


12.3 


1150+497 


Bothwell et al. (2010) 


SMM J131201 


.20+424208. 


.8 


99.60 


0. 


.39 


0.70 


10.2 


1308+326 


Greve et al. (2005) 


SMM J131208 


.82- 


h424129. 


1 


90.62 


0. 


10 


0.60 


18.8 


1308+326 


This work 


SMMJ131232 


.31- 


r423949. 


.5 


103.90 





.39 


2.1 


0.9 


1308+326 


This work 


SMMJ163631 


.47- 


r405546. 


.9 


105.53 


0. 


16 


0.83 


9.7 


3C345 


This work 


SMM J163639 


.01- 


h405635. 


.9 


92.65 


0. 


.28 


1.5 


5.5 


3C345 


This work 


SMM J163650 


.43+405734. 


.5 


102.17 


0. 


07 


0.38 


40.4 


3C345 


Neri et al. (2003) 


SMMJ163655 


.79- 


h405909. 


.5 


95.52 





15 


1.1 


7.0 


3C345 


This work 


SMMJ163658 


.78- 


h405728. 


1 


105.24 





20 


1.0 


8.3 


3C345 


This work 


SMM J163658 


.19- 


H410523. 


.8 


100.11 


0. 


12 


0.73 


14.9 


3C345 


Greve et al. (2005) 


SMM J163706 


.51- 


H405313. 


.8 


102.51 


0. 


.13 


0.76 


13.9 


3C345 


Greve et al. (2005) 


SMMJ221804 


.42- 


h002154. 


.4 


98.32 





.13 


0.63 


24.8 


2223-052 


This work 


SMMJ221735 


.15- 


h001537. 


2 


84.44 


0. 


10 


0.52 


16.5 


2230+114 


Greve et al. (2005) 


SMM J221737 


.39- 


h001025. 


1 


95.52 


0. 


15 


0.78 


14.4 


2223-052 


This work 



RMS noise averaged over the bandwidth of observation. 
RMS noise per channel. 
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Of the 32 SMGs with detected 12 CO emission, it is im- 
portant that we identify any "marginal" cases, for which 
there is a chance that we have a false positive detection. 
For each of our detections, we identify any source with a 
peak flux > 5 times the RMS map noise as a true detection. 
Sources with a peak flux < 5a were individually inspected 
for large spatial and velocity offsets (which we define as > 4" 
from expected phase centre or > 400 km s -1 from expected 
velocity zero point respectively) , either of which result in the 
source being classified as a "candidate" detection. We iden- 
tify four such marginal cases; these are listed in Tables 4 and 
5. In addition, we classify SMM J131208+424I as a 'candi- 
date' detection - despite being just ~ 1" from phase centre 
and ~ 300 km/s from the expected velocity zero point, the 
source is only weakly detected (3.6a), and as such we cannot 
say with confidence that it is a robust detection. These five 
candidate sources are hereafter identified separately in all 
Figures. 

Where we had detectable 12 CO emission, we determined 
the redshift (z) of the 12 CO emission by calculating the 
intensity- weighted peak redshift, 



zco = 



£(*(*)*) 



(2) 



A linewidth can be estimated by measuring the FWHM 
of a Gaussian profile fit to the spectrum. However, if the 
line in question is non-Gaussian (for any number of reasons, 
such as showing an asymmetry, or being double-peaked), 
then a Gaussian profile will be a poor fit. To avoid this 
complication, we use the intensity-weighted second moment 
(s„) of the 12 CO spectrum as a non-parametric estimator of 
the line width: 



J(u — v) 2 I v dv 
J I v dv 



(3) 



where v is the intensity-weighted frequency centroid of the 
line, and I v is the flux as a function of frequency. In most 
cases derived widths are similar (to within the errors) to 
those achieved with traditional Gaussian-fitting. Variations 
occur only in low signal-to-noise galaxies, where Gaussian 
fits struggle to achieve sensible results. 

This non-parametric estimator does have the weakness 
of being sensitive to noise spikes in the spectrum - a large 
spike at the edge of the receiver bandwidth can significantly 
broaden the effective width derived. To reduce this potential 
source of error, we used a Monte Carlo technique to generate 
multiple copies of each spectrum, with Gaussian noise added 
to each (see Table f). We took as our "true" second moment 
the median value of the resultant distribution of intensity- 
weighted second moments. 

The equivalent Gaussian FWHM of the line can then 
be calculated from the second moment, using the relation 



FWHM = 2V21n2 s v ~ 2.35s„ (4) 

The flux density of each observed 12 CO line was found 
by velocity-integrating the 12 CO spectrum, 



^CO = 



I(v)di 



zCO- 2a v 



(for a Gaussian emission line, this is equivalent to inte- 
grating the flux from — 2a to +2a around the centre). The 



12 CO luminosities of the SMGs in our sample were then 
calculated using the standard relation given by Solomon & 



Vanden Bout (2005) 



L^o = 3.25 x 10 7 ScoAv v~ b 2 s Dl (1 + z) 



(6) 



where SAv is the velocity-integrated line flux, Lq Q is the 
line luminosity in Kkms -1 pc 2 , v b s is the observed central 
frequency of the line, and Dl is the luminosity distance in 
Mpc in our adopted cosmology. 

The fits included the option for a positive uniform con- 
tinuum, but only 2 sources appear to have detected con- 
tinuum (see Appendix A and B). The derived limits on 
the continuum emission in the remaining sources are con- 
sistent with those expected from the observed dust spec- 
trum, and exclude significant synchrotron contributions to 
the far-infrared/sub-millimetre continuum in these sources. 

For sources which are not detected in our 12 CO obser- 
vations, we calculate 3-a line intensity limits based on the 
RMS channel noise: 



ICO < 3 RMSchannel \J AV C dv 



(7) 



where RMS c hannoi is the channel noise given in Table 1, 
AVco is the mean linewidth of the detected sample, and 
dv is the bin size in kms -1 . 

Importantly for the interpretation of the evolutionary 
properties of molecular gas, there is no apparent redshift bias 
in the detection probability: a Kolmogorov-Smirnov (K-S) 
analysis of the redshift distributions, sub-mm flux, or 1.4- 
GHz continuum flux suggests that the non-detections are 
consistent with being a random selection drawn from the 
parent sample. 

2.3 Comparing our sample to the wider SMG 
population 

If we wish to use our sample to draw conclusions about the 
SMG population as a whole, it is important to investigate 
whether our 12 CO-detected sample is statistically represen- 
tative of the wider (i.e. 850-/im bright, radio-detected) SMG 
population. For the purpose of comparison, we take as our 
"representative" sample the large spectroscopic database of 
SMGs presented by C05, which functions as an approxi- 
mate parent sample to the sub-mm-selected, radio-detected, 
and optical spectroscopy-confirmed SMGs presented in this 
work. Note that the selection of the C05 sample does imbue 
it with certain unavoidable biases - being sub-mm selected, 
it is biased towards colder sources (see Chapman et al.|2004 



Casey et al.|2009 |, and the radio-selection will tend to bias 



against higher-redshift SMGs with weak radio emission. In 
addition, the necessity of a spectroscopic redshift confir- 
mation will tend to bias the sample against sources which 
are UV-faint, or have weak emission or absorption features. 
Fig. [2] shows a comparison of three observable parameters 
(1.4-GHz continuum flux, 850-/im flux, and J?-band mag- 
nitude), as well as redshift distribution histograms, for the 
detected subset and the C05 parent sample. 

Looking at the distribution of redshifts, our detected 
subset appears to be a representative sample of SMGs (with 
the aforementioned selection effects) at 2 ~ 1-3. We do not 
match the low-redshift coverage, however: low redshift (z < 
1) systems were not included in our programme, and as a 
result the lowest redshift SMG that was targeted in 12 CO is 
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Table 2. List of non-detection observations discarded due to a redshift measurement from the literature lying outside the band, (t SMG 
followed up later at the correct redshift.) 



ID 


z 


z 


Redshift reference 




(targeted) 


(literature) 




SMM J105238.36+572435.3 


3.045 


2.99 


Menendez-Delmestre et al. (2009) 


SMM J131232.34+423949.5t 


2.321 


2.332 


Chapman et al. (2005) 


SMM J123616. 22+621513.2 


2.55 


2.578 


Chapman et al. (2005) 


SMM J123553.33+621337.2 


2.098 


2.17 


Menendez-Delmestre et al. (2009) 


SMM J123635.67+621423.5 


2.005 


2.015 


Swinbank et al. (2004) 



SMM J123629.13 at z = 1.013. Our observations also under- 
sample the higher redshift range z > 3 when compared to 
C05. 

As the remaining three panels show, the observable 
properties of our sample show no significant differences when 
compared to the parent, aside from what results from the 
redshift restrictions. The brightest i?-band sources, for ex- 
ample, lie at z < 1 and do not appear in our sample. The 
850-/im and 1.4-GHz flux distributions of sources appearing 
in our sample are also consistent with the equivalent distri- 
butions for the C05 sample sources, with no apparent bias 
towards more luminous systems. The mean 850^im flux for 
our sample is 6.5 ± 0.6 mjy, compared to 5.9 ± 0.4 mjy for 
the C05 sample (cut to match the redshift coverage - i.e. 
z > 1 — of our sample); the respective radio flux densities 
are 114 ± ll^iJy and 100 ± 12/iJy. KS tests of the 850-^m 
and and 1.4-GHz fluxes distributions suggest that our sam- 
ple is statistically consistent with being drawn at random 
from that of C05. As a result, our sample of 12 CO-observed 
SMGs can be seen as representative of the brighter-end of 
the population of z ~ 1-3 SMGs as a whole, at the era where 
the activity in this population peaks. 



3 CO EXCITATION MODELLING 

In order to derive the gas properties of our SMG sample, it is 
important to fully understand the gas excitation as shown by 
the 12 CO spectral line energy distribution (SLED). The exci- 
tation of a molecular gas reservoir is controlled by the physi- 
cal conditions within the host galaxy: star-forming ULIRGs 
(and SMGs) are expected to have dense, excited gas that 

3 or beyond (e.g. Weiss 



may be thermalised up to J up 

|et al.||2~007| Danielson et al. 2011), in contrast to the low 
density, low-excitation gas which dominates the 12 CO emis 



sion of more quiescent galaxies (e.g. Crosthwaite & Turner 
[2007| pannerbauer et al.||2009| ). 

Even in luminous SMGs, however, the assumption that 
higher- J up transitions always trace fully thermalised gas is 
a poor one - as J up increases, the observed gas becomes 
warmer and denser, and any underlying cold component 



can be missed (see Harris et al. 2010 Ivison et al. 2010 



Ivison et al.|2011 |. If we wish to fully investigate the phys 



ical conditions within our multi-line survey (our data cover 
J up = 2-7), it is important that we fully understand the typ- 
ical SLED of SMGs from low- to high-J up transitions. This 
SLED can then be used to derive a self-consistent excitation 
model. 

Throughout this work, we use O,j-i/io to denote the 



ratio of 12 CO luminositiefj^L' co (J, J - l)/L' co (l - 0) = 
Tf,( J, J — 1) /Tf,(l — 0), where T& is the equivalent Rayleigh- 
Jeans brightness temperature in excess of that of the Cosmic 
Microwave Background. 

In order to directly compare the SLEDs for multiple 
SMGs, it is important to remove the effect of the strong 
dependence of the 12 CO line flux on far-IR luminosity - 
without this step, the variation of L^o with Lfir would in- 
troduce extraneous scatter, and a systematic offset between 
transitions (due to the gradient of the Lq Q -Lfir relation, 
combined with the different luminosities sampled by the dif- 
ferent transitions at different redshifts). We therefore calcu- 
late a "normalising" factor, which will scale the line fluxes 
to the mean Lfir of our sample (= 6.0 ± 0.6 x 10 12 Lq): 

ALco = , T C ° (Lfir— <Lfir>) (8) 
a Lfir 

The gradient of the Lq -Lfir relation has been debated 
in recent years. Efforts at high redshift with small samples 
found a sub-linear slope: 0.62 ±0.08, as given in Greve et al.| 
( 2005 1, for example. Recent analyses have suggested a some- 



what steeper slope, with Genzel et al. (20101 reporting a 



near-linear slope for both SMGs and "normal" star forming 
galaxies. Below (§4.2), we discuss the slope of the Lco~Lfir 
relation for the extrapolated 12 CO (1-0) luminosity, a phys- 
ical quantity that has important implications for the star 
formation efficiency of our SMGs. For the purposes of re- 
moving LFiR-induced scatter, however, it is necessary to ex- 
amine the slope of the relation for the raw (un-extrapolated) 
luminosities of the J up ^ 2 lines. With increasing J up , the 
molecular gas being observed traces progressively warmer, 
denser gas, which is more and more associated with star for- 
mation. The slope of the Lq Q -Lfir relation for the observed 
(i.e. J up ^ 2) 12 CO luminosities will therefore be dependent 
on Jup. 

Figure |3] shows the Lc -Lfir correlation for our sam- 
ple of SMGs, colour-coded by the observed transition. We 
find that the slope does indeed steepen with increasing 
Jup, finding slopes of dL' co ^ 2 -i) /dLFm = 0.71 ± 0.19, 
^Lco(3-2)/^Lfir = 1-2 ± 0.3, and dL CO ( 4 _ 3 )/<f Lfir = 
2.1 ± 0.4. Of course there is the potential for a redshift- 
induced bias in this result - the J up = 2 lines are observed 
to lower redshifts than the more excited transitions, which 
could result in the artificial "flattening" of the relation. We 
do note, however, that the resulting line luminosity ratios 
reached are relatively insensitive to the choice of slope - 
adopting a shallower slope has the effect of lowering the line 
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Figure 2. Plots comparing the observed properties of the 1 CO sample presented in this work, along with the "parent" sample of C05. 
Here (and in plots hereafter), open circles denote the SMGs we identify as "candidate" detections, while filled circles are our secure 
detections. All four panels have redshift as the abscissa. The panels: Top left: 850-/xm flux. Top right: 1.4-GHz continuum flux. Bottom 
left: ij-band magnitude. Bottom right: Redshift distributions. Sources comprising our sample do not differ significantly from those in 
Chapman et al. (2005), over a similar redshift range. 



luminosity ratios, but for a reasonable range of slopes (i.e. 
> 0.6) the variation is within the bootstrap errors. 

After the removal of this Lpir dependence, SMGs with 
Lfir = < Lfir > will have their line fluxes unchanged, and 
SMGs which are less (more) luminous will have their line 
fluxes proportionally increased (decreased). 

The resulting 12 CO-line SLED for our sample of SMGs, 
adjusted to a mean far-IR luminosity, is shown in Fig. [4] 
The fluxes have also been adjusted to a common source red- 
shift of z = 2.2. In addition to the 12 CO line observations 
presented in this work, we have compiled all the observed 
transitions available in the literature for SMGs in our sam- 
ple (compiled from |Carilli et al.|201 Ricchers et al. (20101; 
Ivison et al. 2011). SMGs with multiple- J up transitions are 
highlighted in Fig. [4] and we also show the trend of the me- 
dian 12 CO luminosity/flux as a function of J up , which rep- 
resents the SLED of "typical" SMGs. 

As can be seen from the median line flux in the top 
panel of Fig. [4] the 12 CO SLED shows a moderate excitation. 
There is a monotonic increase in Sco up to J up ~ 5, with 
a flattening or turn-over above this. This is similar to local 



star burst galaxies (see, e.g. Weiss et al. 20071. The behaviour 



of the best sampled SMG in our sample (SMM J123711.86 
at z = 4.04 has five observed transitions) is not dissimilar to 
that of the median, showing an increase up to the J up = 5 



transition, and a turnover at higher excitations. We do cau- 
tion, however, that the behaviour of the median SLED above 
the apparent turnover is poorly sampled, and dominated by 
low number statistics; hence constraints on the higher exci- 
tation components must be taken as tentative. 

The middle panel of Fig. [4] shows our derived 12 CO 
SLED, compared to three other well-studied systems, nor- 
malised to the 12 CO (3-2) flux. The moderate excitation 
found for our SMGs is certainly more excited than less- 
active, local galaxies. It is very similar, however, to the well- 
studied, stron gly lensed SMG SMM J2135— 0102 ("The Cos- 
mic Eyelash" ; Danielson et al.|2011 1 , showing a similar peak 
and falloff towards higher- J up . 

The bottom panel of Fig. [4] shows the results of fit- 
ting our SLED with the Photon Domimated Region (PDR) 
model of Meijerink et al. ( 2007[ ). As can be seen, a single- 
component model provides a very poor fit to the data. Hence 
we are driven to a two-component model, consisting of sep- 
arate "warm" and "cool" components, which provides a 
much better fit to the data. The best fit is achieved with 
a combination is with 75% of a cool component plus 25% 
of the warm component. The cool component has a den- 
sity log 10 (n) = 2.0cm -3 and log 10 (G) = 1.5 (where n is 
the Hydrogen volume density density, and G is the far-UV 
field strength in Habing units), and the warm component 
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Figure 3. The Lqq - ^fir relation for our sample of SMGs, 
colour-coded by the observed transition. Slopes are fitted to the 
respective samples. We find that a fit to the J up = 2 points has a 
slope of 0.7f ± 0.19, the J up = 3 points have a slope of 1.2 ± 0.3, 
and the J up = 4 points to have a slope of 2.1 ± 0.4. Hence the 
LJ-jo - Lfir relation is progressively steeper for increasing values 
of Jup, as the observed gas becomes a more efficient tracer of 
active star formation. 



has a density log 10 (n) = 5.5cm -3 and a radiation field of 
log 10 (G) = 2.0. Thus the volume-averaged radiation field 
experienced by both components is between 30-100 x that 
of the Milky Way. Given the high star-formation rates for 
our SMG sample (~ 1000 x the Milky Way), these relatively 
modest radiation fields suggest their star formation must be 
spatially extended. Similarly, the characteristic densities we 
derive range between that expected for the ISM in typical 
star-forming galaxies and that in dense starburst systems 
(see Fig. 4 in Danielson et al. 2011). As the bulk of the mass 
in these galaxies is in the cool component, we can use the 
best-fit density for this component and assume the gas is in a 
1-kpc thick disk with a radius of R ~ 3 kpc (see §4.4) to pre- 
dict a typical gas mass for an SMG of the order of ~ 10 11 M 
comparable to the masses we estimate later in §4.4. However, 
we caution that acceptable model fits to the SLED span a 
significant range in parameters: the cool component can have 
a density of 2.0 < log 10 (n) < 3.0 and a range of acceptable 
radiation fields of 1.0 < log 10 (G) < 2.5, while the warm 
component parameters have ranges 4.5 < log 10 (n) < 6.5 
and 2.0 < log 10 (G) < 3.0. 

Using the SLED shown in Fig. [4] we calculate median 
brightness temperature ratios (equivalent to line luminosity 
ratios) which we use to convert our J up ^ 2 observations into 
an equivalent 12 CO J = 1-0 flux. These are given in Table[3j 
and we use them throughout this work. These values agree 
well with values reported elsewhere in the literature (e.g. 
Ivison et al. 2011). 



4 THE PHYSICAL PROPERTIES OF SMGS 

In this section, we describe the modes and methods used to 
derive physical parameters for our sample of SMGs - these 
are given in Table 5. 
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Figure 4. The 12 CO SLED for the SMGs in our sample. All 
SMGs have been normalised to the mean far-IR luminosity of 
the sample as detailed in the text, and fluxes have been adjusted 
to the median redshift of the sample (z = 2.2). Top: The 12 CO 
SLED for our SMGs, in (normalised) flux units. The SMGs in 
our sample with multiple- J up observations available in the liter- 
ature are shown by connecting dotted lines. The colour-coding 
is shown in the inset legend. SMGs with only a single observed 
transition are shown as blue circles. The median SLED of the 
sample, along with the bootstrapped error, is overlaid. Middle: 
The median SLED compared to other well-studied galaxies, nor- 
malised to the flux in the 12 CO (3—2) transition. Bottom: The 
median SLED, fit with PDR models as discussed in the text. The 
median SLED shows a moderate excitation, similar to that of 
SMM J2135-0102 and is best fit by a two-component model of 
the ISM. 
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Table 3. Median brightness temperature ratios for the SMGs in 
our sample. All errors estimated using bootstrap resampling. 



Transition 



Value 



'"21/10 
r 32/10 
r 43/10 
r 54/10 
r 65/10 
r 76/10 



0.84 ± 0.13 
0.52 ±0.09 
0.41 ± 0.07 
0.32 ±0.05 
0.21 ± 0.04 
0.18 ±0.04 
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Figure 5. Distribution of Gaussian-equivalent FWHMs for the 
sourced detected in 12 CO. FWHMs arc derived from the intensity- 
weighted second moment, as described in §4.2| SMGs detected in 
"high-J U p" (J up ^ 3) and "low-J up " (J up = 2) are plotted. As 
we discuss in the main text, we believe that the lower FWHM 
displayed by the lower-J up 12 CO lines are likely to be due to 
their lower mean redshifts and hence lower luminosities. 



4.1 



"CO Kinematics 



The 12 CO line emission from SMGs traces the kinematics 
of the potential well in which the molecular gas lies. Previ- 
ous studies of this emission have uncovered very broad lines 
Greve et aL]|2005| Tacconi et al. 2006, 2008): |Greve et al 



2005) found a mean FWHM linewidth of 780kms _1 for 



their SMG sample, suggestive of deep gravitational poten- 
tial wells. 

The mean FWHM of our sample (which includes the 
Greve et al. 2005 SMGs) is 510±80kms~ 1 . This is somewhat 



lower than that found by Greve et al. ( 2005 1 and is most 



likely attributable to the bias towards IR-luminous sources 
in Greve et al. ( 2005 1 - the far-IR luminosities of their sam- 
ple are higher than our sample, by a factor of ~ 2.5 (they 
quote < Lpir >= (15 ± 7) x 1O 12 L0, whereas our sample 
has <L F tr>= (6.0 ±0.6) x 10 12 L Q ). As there is a strongly 
positive Lfir -Mdyn correlation (see Fig. |6j , it is to be ex- 
pected that the more IR-luminous SMGs presented in | Greve] 
et al. ( 2005| ) would have broader line widths. 

The distribution of FWHM values for our sample is 
shown in Fig. [5] The values derived from low- J up (i.e. J = 2- 



1) and higher- J up (i.e. J up ^ 3) lines have been highlighted 
separately. There is a slight bias towards lower line widths 
in the low- J up subset: the J = 2-1 observations have a mean 
FWHM of 470 ± 80 km s~\ compared to 550 ± 90 km s" 1 
for the J up ^ 3 sample. This trend is the opposite of what 
might be expected (everything else being equal, the lower- 
Jup lines trace more extended gas than the higher- J up lines, 
and therefore the linewidth should reflect a higher dynami- 
cal mass), and suggests a redshift-induced bias in the pop- 
ulation, with low-Jup transitions being measured in lower- 
redshift, typically lower-luminosity SMGs. 

Figure |6j left panel, shows the FWHM of the observed 
emission line plotted against the (radio-derived) far-IR lumi- 
nosity. Also shown on the plot are the "submillimetre-faint 



radio galaxies" (Chapman et al. 2008; Casey et al. 20111, 
warmer dust counterparts of SMGs, and z ~ 2 optically se- 
lected star-forming galaxies (Tacconi et al. 2010), less active 
star forming galaxies at these epochs. 

Of the three samples of galaxies shown in Fig. [6j our 
SMGs clearly have the largest linewidths. There is a signifi- 
cant overlap between the SFRGs and the SMGs, suggesting 
that the "hotter dust" ULIRGs are kinematically similar 
to SMGs. The less-actively star-forming galaxies (SFGs), 
however, have only a slight overlap with the SMGs, hav- 
ing typical linewidths < 300 kms -1 - lower than all but 
the most narrow-line SMGs. Interestingly, the SFGs have 
far-IR luminosities a factor of several lower than even those 
narrow- line SMGs/SFRGs with compatible 12 CO FWHMs, 
suggesting that even SMGs which are kinematically com- 
parable to optically-selected galaxies have enhanced star 
formation. Dissecting the mechanisms driving this strong 
emission in the narrow-line SMG population requires high 
resolution 12 CO imaging, in order to spatially resolve the 
kinematics. 

Having measured 12 CO line luminosities and line widths 
for our sample, we now turn to the correlation between these 
two observables, which respectively trace the mass of the gas 
reservoir and the dynamics of the potential well in which it 
lies. Figure [©J right panel, shows the FWHM of the 12 CO 
emission line plotted against the derived 12 CO (1-0) lumi- 
nosity for our sample of SMGs. Also shown on the plot 
are the 12 CO (1-0) observations for SMGs by Ivison et al. 
(2011). Included on the plot (but not included in any of 
the fits) are the z ~ U/LIRGs from Downes & Solomon 
(1998). We fit a power-law fit to the SMG sample, deriving 
a power-law index of 1.98 ± 0.07. The scatter around this 
power-law is low - just AL' co /L' co = 0.38. 

As L' co is sensitive to the total mass of molecular gas, 
while the FWHM is sensitive to both dynamical mass and 
any inclination effects, some dispersion around the relation 
would be expected. A population of thin, randomly orien- 
tated discs would introduce an inclination-based dispersion 
of ~ 50% - the fact that our overall scatter is lower than this 
expected value strongly suggests that we are instead seeing 
gas in the form of thick discs (or turbulent ellipsoids), for 
which the line of sight velocity depends less strongly on in- 
clination. 

We also overlay a simple functional form for Lq O ( 1 _ j: 



C(V/2.35) 2 i? 



(9) 



CO(1 - 0) 1.36aG 
where V is the FWHM 12 CO line width; 1.36a is the conver- 
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Figure 6. Left panel: LpiR plotted against FWHM for the sample of 12 CO-detected SMGs. Also plotted are SFRGs from Chapman et 
al. (2008) and Casey et al. (2011), and z ~ 2 SFGs from Tacconi et al. (2010). The SFGs have IR luminosities several times lower than 
SMGs and SFRGs with compatible linewidths. The blue dashed line shows a power-law fit to the SMGs alone (index = 2.2 ± 0.2), and 
the black dotted line shows a fit to all three samples of galaxies (index = 2.6 ± 0.3). The SMGs have broader typical linewidths than 
other populations, and have enhanced far-IR luminosities compared to SFGs (even at comparable linewidths). Right panel: Lco(i— 0) 
plotted against FWHM for the same sample. Also included are values for local U/LIRGs from Downes & Solomon (1998), and J =1—0 
observations of SMGs from the literature, but these are not included in the fits. The solid purple line shows the relation described in 



f 4.1 



while the dashed black line shows a simple fit to the points. The SMGs show a tight correlation between Lco(l-O) an d FWHM. 



sion factor from 12 CO line luminosity to gas mass (we adopt 
a = 1; see §4.4 below), C is a constant parameterising the 
kinematics of the galaxy (where we have taken C — 2.1 ap- 
propriate for a disc: see Erb et aL"||2006 ), 7? is the radius of 
the 12 CO J = 1-0 emission region which, following Ivison et 



al. (2011) we have taken to be 7kpc (see f 4.2 below for dis- 
cussion of the issue of 12 CO sizes), and G is the gravitational 
constant. 

This simple model provides a very good fit to the ob- 
served trend. Adopting Lq O( . 1 _ j oc V 2 is clearly a good 
match to luminous SMGs. The tight correlation seen in high- 
redshift ULIRGs is in stark contrast to z ~ U/LIRGs, 
which show very little correlation between V and Lq O j 1 _ j. 
We suggest that this is due to a combination of different ge- 
ometry of the gas reservoirs; the thin nuclear gas discs and 
rings seen in local U/LIRGs mean the lack of any inclination 
correction produces significant scatter in our plot, combined 
with a range in / gas and a more significant stellar mass con- 
tribution to the dynamics of the potential well (there is little 
or no correlation between the stellar masses of the SMGs in 
our sample and their 12 CO linewidths). 

The strength of the Lq -FWHM correlation either in- 
dicates a very uniform ratio of gas-to-stellar contribution to 
the dynamics of the region probed by the 12 CO emission, or 
that the gas mass dominates in this region (although that 
then requires the stellar mass to be significantly extended 
beyond the gas radius probed with our observations). There 
is little evidence for this, however; Swinbank et al. (20101 



report Hubble Space Telescope NICMOS ff-band radii of 
2.8 ±0.4 kpc for the SMGs, comparable to the extent of the 



4-1.1 Double-peaked sources 

Looking in more detail at the spectra, some SMGs exhibit 
double-peaked 12 CO spectra, a potential indication of kine- 
matically d istinct components within these systems. |Greve| 
et al.| < |2005| > found that 4-6 of their 12 12 CO-detected SMGs 
displayed evidence of multiple kinematic components (as did 
Tacconi et al. 2006, 2008 and Engel et al. 2010), and also 
suggested these tended to be the lower-redshift systems. Our 
larger sample allows us to investigate these trends in more 
detail: We find that 7 of the 32 12 CO-detected SMGs in the 
full sample are clearly double-peaked, and a further 2 have 
spectra that are better-fit by a double than a single Gaus- 
sian profile (measured by the reduced x 2 of the fits). This 
is a lower fraction, 20-28%, than previous findings. In ad- 
dition we do not find any evidence that the double-peaked 
sources lie at systematically lower redshift than the sam- 
ple as a whole, in contrast to the suggestion of Gre ve et al.| 
(20051. 



high- J U p 12 CO emission jEngel et aL]|2010| |Bothwell et al 



2010) 



It must be noted that our observations can fail to de- 
tect kinematically-distinct components in two situations. 
Firstly, if the velocity separation between the compo- 
nents is too small they will appear as a single line. The 
typical peak separation of the double-peaked SMGs is 
~470kms _1 , comparable to the mean line width of the 
sample (~ 500 kms -1 ) - velocity separations much smaller 
than the line-width will be blended. Secondly, the low spa- 
tial resolution of our maps cannot distinguish spatially- 
separated components with similar line-of-sight velocities; 
our source SMM J094303+4700, for example, has been ob- 
served with higher-resolution imaging, finding two distinct 
components (termed H6 and H7; Ledlow et al. 2002). This 
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is also true of SMM J123707+6214 (Tacconi et al. 2008) and 
SMM J 105307+5724 (Bothwell et al. 2010). 

Thus our 20-28% fraction of double-peaked profiles 
should be interpreted as a lower limit to the true rate in 
SMGs. 

Examining the physical properties of the double-peaked 
sources, we find that they do not seem to differ significantly 
from the population as a whole. The mean gas mass for 
double-peaked spectra is (6.2 ± 1.2) x 10 10 Mq, compared 
to (5.3±1.0) x 10 10 Mq for the complete sample. The mean 
far-IR luminosities for the respective classes differ by a com- 
parable about, with the double-peaked SMGs having a mean 
Lfir = (7.2± 1.1) x 10 12 Lq, while the single-peaked galax- 
ies have a mean Lfir = (6.0 ± 0.6) x 10 12 Lq. That is, 
double-peaked sources have molecular gas reservoirs and far- 
IR luminosities approximately 20% times greater than those 
with single peaks, but given the small numbers of sources 
examined these results are certainly not significant. 



4.2 Dynamical Masses 

Our measurement of the kinematics of the 12 CO emission 
from SMGs allows us to estimate masses for the galaxies 
themselves from the width of the 12 CO line, with an as- 
sumption about the dynamical structure and extent of the 
system. For example if the 12 CO emission arises in a vit- 
alised body of radius R, with one dimensional velocity dis- 
persion a = s u , then the dynamical mass is given by: 



Mdyn [Mq] 



(10) 



An alternative model for calculating dynamical masses 
is to assume that the line emission originates from a rotating 



disc. In this case, the dynamical mass is given by (i.e. Neri 
et al.|2003| |: 



M 



dyn • 



[Mq] 



4 x 10 4 V 2 R, 



(11) 



where V is the FWHM velocity, and i is the inclination 
of the disc. For the purpose of calculating dynamical masses 
using this method, we adopt a mean inclination angle, ap- 
propriate for a population of randomly-orientated discs, of 
< sini >= 7r/4 ~ 0.79 (see Appendix A of Law et al.|2009 |. 

The choice of which estimator to use depends on the as- 
sumed form of the molecular ISM. If we assume that bright 
SMGs are predominantly major merger events (as argued 
by Engel et al.||2010 l then the virial estimator may be the 
correct choice. However, recently the dynamics of the 12 CO 
emission from some SMGs (such as the lensed z = 2.3 SMG, 
SMM J2135— 0102, which has been mapped at 100-pc reso- 
lution; Swinbank et al. 2010, 2011), have revealed that the 
molecular gas appears to lie in a rotating disc. However, the 
intrinsic luminosity is at the fainter end of our sample (Ivi- 
son et al. 2010). This is consistent with the suggestion made 
above that there are a mix of disc-like and virialised systems 
in the SMG population, with discs being more common at 
the lowest luminosities, while the brighter end of the SMG 
population has a higher frequency of merging systems. Of 
course, a rotating disc configuration does not preclude the 
conclusion that the galaxy is undergoing a merger, though it 
does suggest that any such merger is likely to be late-stage. 

The final source of uncertainty in these calculations 
is the extent of the 12 CO emission. While the compact- 
configuration PdBI observations in our programme are well- 



suited to a detection survey, they have the disadvantage of 
lacking the angular resolution necessary to resolve the emis- 
sion in typical SMGs. Higher-resolution observations from 
PdBI (typically in the higher- J up transitions) have suggested 
that the warm gas reservoirs in SMGs have measured radii 
of 2-3 kpc jTacconi et~aT1|2008||Bothwell et al.||2010| |Engel 



et al. 20101. Ivison et al. (2011), using ^CO(l-O) obser- 



vations of SMGs, suggested that the derived radius is also 
dependent on the transition observed; higher- J up lines pref- 
erentially trace denser and more centrally concentrated star- 
forming gas, yielding smaller radii, than the more extended 
reservoirs traced by lower- J up emission (which includes a 
more-extended component which may not be directly asso- 
ciated with vigorous star formation). 

Examining all SMGs observed in high-resolution 12 CO 
(see the compilation presented by |Engel et al.|[2010| ), there 
seems to be something of a transition-based effect on the 
observed size of the 12 CO reservoirs in SMGs; the galaxies 
observed in 12 CO (4-3) emission have a mean radius of 2.1 ± 
0.7 kpc, while those observed in 12 CO (3-2) have a mean 
radius of 3.1 ± 1.4 kpc (though these results should be taken 
as somewhat tentative, as sample sizes are small here). The 
redshift ranges covered by the respective samples are similar, 
suggesting that this is not an evolutionary effect. These high- 
J up radii are smaller than the extents measured for 12 CO (1- 
0) emission; Ivison et al. (2011) found typical radii of ~ 7 kpc 
for a smaller sample of five SMGs, observed with the JVLA. 
Those authors also find some evidence that the FWHM of 
the 12 CO (1-0) emission is ~ 15% broader than the higher- 
J up lines measured for the same SMGs, roughly consistent 
with the size variations. 

With a likely variation in radius with transition, it could 
be argued that the correct radius to use for the dynamical 
calculation is the one appropriate for the transition used to 
derive the FWHM. However, one of the reasons for deriving 
the dynamical masses is to compare them to the "total" gas 
and baryonic masses of these galaxies, we need to bear in 
mind that these "total" measurements are effectively mea- 
sured in an aperture equal to the size of the 12 CO (1-0) 
emission (as we are using galaxy-integrated O.j-i/io val- 
ues). The JVLA results suggesting a larger 12 CO radius are 
however, as yet, based on small sample sizes. 

To avoid introducting more uncertainty due to the 
choice of assumed size (which our low-resolution survey is 
not optimised to study), we explicitly retain the size depen- 
dence in our estimates of the dynamical masses. (For ease of 
comparison with other studies, we shall also quote the final 
values assuming both "extreme" cases) 

The median dynamical mass derived for our sample 
depends on the choice of mass estimator. Adopting the 
"virial" estimator given in Eq. 6, the median dynamical 
mass of the sample is (7.1 ± 1.0) x 10 10 R kpc MqF] Re- 
moving the dependency on radius, this value corresponds 
to (2.1 ± 0.3) x 10 11 Mq if we assume R = 3 kpc, and 
(5.0 ± 0.7) x 10 11 Mq if we assume R = 7 kpc. 

Alternately, adopting the "rotational" (i.e. disc-like) es- 
timator, given in Eq. 7, results in a median dynamical mass 
of (1.6 ± 0.3) x 10 10 R kpc M . A gain removing the depen- 



3 The error here is the statistical error resulting from uncertainty 
on the FWHM measures. 
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Figure 7. The star formation efficiency, LJ-, Q vs. Lfir. fncluded 
on the plot are two samples of local ULIRGs, and the J=(l — 0) 
SMG observations of Ivison et al. (2011). Best fitting slopes to 
the SMGs alone (blue dotted line), the local ULIRGs alone (grey 
dot-dashed line), and all three combined samples (red dashed line) 
arc ovcrplottcd. They have slopes of 0.93 ± 0.14, 0.79 ± 0.08, and 
0.83 ± 0.09 respectively. 



dency on radius, this value corresponds to (4.8 ± 0.9) x 
10 10 M if we assume R = 3 kpc, and (1.1±0.2) x 10 11 M 
if we assume R = 7 kpc. 

Each of these results are roughly in line with what might 



be expected from total halo masses; Hickox et al. (20121 
found, using a clustering analysis, that SMGs typically lie 
within Dark Matter haloes of masses Mhaio = 9 x 10 12 Mq. 

These dynamical masses are also reasonably consistent 
with recent estimates for SMGs derived using other indica- 
tors. Swinbank et al. (2004; 2006) used the kinematics of 
the spatially-resolved Ha emission line to estimate the dy- 
namical mass of a sample of bright SMGs, finding that a 
mass of (5 ± 3) x 10 11 Mq represented the ensemble pop- 
ulation well. The high dynamical masses found for SMGs 
are, in general far greater than those found in other, less 
extreme samples of galaxies: galaxies comprising the UV- 
selected SINS sample, for example have a typical dynamical 
mass of 8 x 10 10 M ( Forster Schreiber et al.|2006||Shapiro| 
et al.|[2009 l. Due to uncertain model parameters (size, and 
choice of kinematic structure) inherent to a low-resolution 
study such as this, the range of possible dynamical masses is 
large. It is worth noting that a comparison to values derived 
from other studies favours the upper end of our estimates, 
suggesting that systems which are both compact and disc- 
like perhaps do not comprise a significant proportion of the 
SMG population, and that extended and/or merging SMGs 
represent most of the SMGs in current studies. See §4.6| for 
more discussion of this issue. 



4.3 



The Lq — Lfir correlation 



A useful quantity to measure for our sample of SMGs is the 
efficiency with which their molecular gas is being converted 
into stars. The star formation efficiency (SFE) is sometimes 
defined as SFR/M(H2) - the inverse of the gas depletion 



time - but here we take the approach of parameterising this 
as a ratio of observable quantities: Lfir and Lco(i-o)- 

There has been debate, in recent years, as to the value of 
the slope of the Lco(i-o) - Lfir relation. Whereas earlier we 
discussed the difference in slopes between the various tran- 
sitions observed (in order to derive the median SLED), here 
we present the relation between the derived 12 CO J =1-0 
luminosity and Lfir - a relation which describes, in observ- 
able terms, the relationship between the luminosity due to 
star formation, and the total gas content. 

In their study of 12 SMGs, |Greve et a\] pOOl] ) found 
a slope of the relation between Lfir and L c 



CO(1-0) 



of 



0.62 ± 0.08 fit a combined sample of lower-redshift LIRGs, 
ULIRGs, and SMGs - identical to the slope derived for the 
local LIRGs/ULIRGs alone. The small number of galaxies in 
Greve et al. (20051, however, prevented a full investigation 



of the SFE slopes within the SMG population itself. Some 
recent authors, however, have found the slope to be closer 
to linear -|Genzel et~aL] (|2010|> found a slope of 0.87 ± 0.09 



fit a combined sample of star-forming galaxies across a wide 
range of redshifts. 

Figure [7] shows the Lq Q — Lfir relation for our sam- 
ple of SMGs. Included on the plot are datapoints for lo- 



cal (U)LIRGs, as measured by Sanders et al. ( 1991 1 and 
| Solomon et al.| ( |1997[ ). We also show three power-law fits; 
to the local (U)LIRGs alone, the SMGs alone, and the 
combined sample. We find the SMGs to lie slightly above 
the best fitting (sub-linear) line for local (U)LIRGs, ne- 
cessitating a steeper slope. The power-law fit to the local 
(U)LIRGs alone has a slope of 0.79 ± 0.08, while the fit to 
the combined sample of SMGs and (U)LIRGs has a slope of 
0.83 ± 0.09 - very close to linear. It can also be seen that 
a fit to the SMG sample alone has an even steeper slope 
of 0.93 ± 0.14 - although, within the uncertainty, this is 
consistent with the slope for the combined samples. These 
results are in good agreement with most previous findings; 
the near-linear slopes (which agree well with those found 
by Genzel et al. 2010) would imply a roughly constant gas 
depletion timescale across the entire range of far-IR lumi- 
nosities shown here. 

However, we caution that our analysis requires extrap- 
olating from high-Jup 12 CO transitions to 12 CO(1-0). As 
we showed in Fig. 3, plotting Lc -Lfir for different tran- 
sitions shows increasingly steep slopes with higher J up and 
hence using a single brightness temperature ratio to cor- 
rect each of these correlations to a "common" transition 
will introduce significant scatter and uncertainty into the 
resulting slope. Indeed, Ivison et al. (2010) have under- 
taken a similar analysis based solely on 12 CO (1-0) obser- 
vations and homogeneously-derived far-IR luminosities and 
conclude that the Lq O ^ 1 _ j-Lfir has a slope substantially 
below unity (consistent with the trend of slope with J up seen 
in Fig. 3). We therefore suggest that it is difficult to draw 
any strong conclusions from high- J up observations about the 
gas depletion timescales of the different populations or the 
form of the Kennicutt- Schmidt relation in these galaxies, 
as these are too uncertain without brightness temperature 
ratio measurements for individual sources. 
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4.4 Molecular gas masses 

Part of the power of observations of 12 CO emission from 
high-redshift galaxies is that they provide a powerful tool 
to derive the mass of the reservoir of molecular gas in these 
systems, which is mostly in the form of H2. This is of critical 
importance because this reservoir is the raw material from 
which the future stellar mass in these systems is formed. 
Along with the existing stellar population, it therefore gives 
some indication of the potential stellar mass of the resulting 
galaxy at the end of the starburst phase (subject, of course, 
to the unknown contribution from in-falling and out-flowing 
material). 

Estimating the mass of H2 from the measured L^o re- 
quires two steps. Firstly, luminosities originating from higher 
transitions (J up ^ 2) must be transformed to an equivalent 
12 CO J = 1-0 luminosity, using a brightness ratio. We have 
derived the necessary brightness ratios using our composite 
SLED as discussed in [J3] above. Once a Lq O ( 1 _ j has been 
determined, it must be converted into a H2 mass by adopting 
a conversion factor a: M(H 2 ) = aL' co , where a is in units of 
Mq (Kkms" 1 pc 2 ) -1 (when discussing a hereafter, we omit 
these units for the sake of brevity). This can then be con- 
verted to a total gas mass, including He, M gaa = 1.36M(H2). 

There is a large body of work, both observational and 
theoretical, dedicated to determining the value - and ascer- 
taining the metallicity or environmental dependence - of a 



(e.g. Young & Scoville|1991 


Solomon & Vanden Bout 


2005 


Liszt et al.||2010| |Narayanan et al.|2011| |Genzel et al. 


2011 


Papadopoulos et al.||2012 


I. While secular discs such as the 



Milky Way have a relatively "high" value of a ~ 3-5, using 
this value for the gas in nuclear discs/rings within merging 
systems and starbursts at z ~ leads to the molecular gas 
mass sometimes exceeding their dynamical masses. As such, 
a lower value - motivated by a radiative transfer model of the 
12 CO kinematics - is typically used for the intense nuclear 
starbursts in the most IR-luminous local systems: a ~ 0.8, 
with a range of 0.3-1.3 (Downes & Solomon 1998). However, 
some recent results have suggested that this value might, in 
fact, underestimate the true value in high-redshift SMGs. 
Bothwell et al. (2010) found that applying the canonical 
ULIRG value to two 2 ~ 2 SMGs resulted in gas fractions 
of < 10%, incongruous with their extreme star formation 
rates. Similarly, a dynamical analysis has been undertaken 
on the high-redshift SMG, SMM J2135-0102, by Swinbank 
et al. (2011) yielding a higher value, a ~ 2 (supported by 
LVG modelling; Danielson et al. 2011). 

Here we adopt a value of a =1.0, and caution that all 
gas masses derived are dependent on this uncertain param- 
eter. Using this value, the resulting mean H2 mass of our 
sample SMGs (including limits from the non-detections) is 
(5.3 ± 1.0) x 10 10 Mq. This is comparable to the findings 



of |Greve et al. (20051, once the conversion factor and exci- 
tation model (detailed in §j3j> is accounted for; adjusted for 
our model values, the SMGs presented in that work have a 
mean gas mass of (6.5 ±2.2) x 10 10 Mq. 

The distribution of H2 masses for our sample is shown 
in Fig. [8] As previously, we have differentiated in the distri- 
bution between low-J up (J up = 2) and higher- J up (J up > 3) 
lines. There is a notable tendency for the low- J up lines to re- 
sult in lower gas masses than the higher- J up subsample. The 
mean gas mass for the J up = 2 subsample is (3.2±2.1) x 10 10 
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Figure 8. Distribution of H2 gas masses for the full sample. To 
calculate H2 masses, a 12 CO/H2 conversion factor of a =1.0 has 
been assumed as detailed in j ]4.4| We identify both "high-J up " and 
"low-J up " derived masses; the offset between these most probably 
arises due to a redshift bias, whereas the "low-J up " sources are 
smaller galaxies lying at a lower mean redshift. 



Mq , compared to a mean gas mass of (6.2 ±0.2) x 1O 1O M 
for the SMGs observed in higher J up transitions. This factor 
of ~ 2 difference is substantially more than could be at- 
tributed to random errors in the adopted brightness model, 
and it likely a result of the low- J up observed SMGs lying at 
a lower mean redshift (and therefore having a lower mean 
far-IR luminosity; see Fig. [5] where a similar effect is seen 
in the distribution of linewidths). Splitting our SMG sam- 
ple into low-redshift (2 < 2) and high-redshift (2 > 2) sub- 
samples, we find H2 masses of (4.1 ±1.0) x 10 10 Mq and 
(6.1 ± 0.2) x 10 10 Mq respectively. 

The gas masses determined for our sample of SMGs are 
roughly comparible to those derived for other high-redshift 
star forming galaxies (though differences in 12 CO/H2 con- 



version factor make comparisons difficult). Tacconi et al. 
(20101 find mean molecular hydrogen masses of (7.9 ±2.9) x 



10 1U M and (1.3 ± 0.7) x 10 11 M respectively for their 
low-redshift (2 ~ 1.2) and high-redshift (2 ~ 2.3) samples 
of "main sequence" star-forming galaxies (selected from the 
AEGIS and BX/BM surveys). These are similar to our es- 
timates, but we note that the two samples were calculated 
with very different values of a - Tacconi et al. (20101 adopt 
a value of a = 3.2, more appropriate for secular disc galaxies 
and a factor of three higher than the a = 1 we use for our 
SMGs. Using a similar choice of a for both samples would 
obviously lead to disagreements in their gas masses. 



4.5 The evolution of the gas fraction 

If it is true that all star formation follows a universal scaling 
law dependent solely on the available gas, then the peak in 
cosmic star-formation activity - and subsequent decline - is 
simply a reflection of the availability of molecular gas within 
galaxies. Recent results have strengthened the expectation 
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Figure 9. The redshift evolution of the baryonic gas fraction. 
The evolution of the mean baryonic gas fraction for our SMG 
sample, calculated in bins of width Az = 0.5, is overlaid (due 
to the paucity of high-redshift SMGs the final bin is the mean 
for all z > 3 SMGs). We also plot star-forming galaxies from the 
literature as follows: at z ~ (Bothwell et al. 2009); z ~ 0.4 
| |Geach et al^|2l)il) ; z ~ 1.5 (Daddi et al. 2010); z ~ 1.2 and 
z ~ 2.3 (Tacconi et al. 2010). We also plot, as a grey shaded area, 
the GALFORM predictions for the evolution of the baryonic gas 
fraction for LIRGs (LprR. > 10 Ln) with halo masses between 
10 12 and 10 13 Mq . There appears to be very little evolution in 
the mean SMG gas fraction over a wide interval of cosmic time. 
Furthermore, our SMGs exhibit similar gas fractions to lower lu- 
minosity, optically-selected populations at comparable redshifts. 
A schematic representation of the evolution of the star formation 
rate density (Sobral et al. 2012) is overlaid - it can be seen that 
the behaviour of the SFRD broadly follows the global evolution 
of the gas fractions. 



that normal, star-forming galaxies at high redshift are sub- 
stantially more gas rich than their z ~ counterparts, by 
perhaps a factor of 3-10 ( Tacconi et al.pOlO Daddi et al 



2010 Geach et al. 2011). As pointed out by Dave et al 
( |2011[ ), the drop in gas fraction between high-redshift and 
the present day could be a straightforward result of the gas 
supply rate dropping faster than the gas consumption rate. 
Conversely, a gas fraction which rises with time represents 
a rapid accretion rate which swamps the consumption rate, 
leading to the accumulation of large gas reservoirs (thought 
to occur at z ^ 4; |Papovich et al.pjll) |Dave et al.||2011[ ). 
The implication of this, of course, is that the enhanced star- 
formation rates seen in these normal star-forming galaxies 
at z > 1 do not necessarily represent a "super-efficient" 
mode of star formation, but could instead simply reflect the 
larger gas reservoirs available in the early Universe (poten- 
tially driven by correspondingly larger accretion rates from 
their surroundings). 

The gas content of galaxies is controlled by three main 
processes; gas accretion via inflow from the intergalactic 
medium (IGM), gas consumption (and associated outflows) 
driven by star formation (and perhaps AGN), and out- 
flows. We can define the baryonic gas fraction, / gas — 
M gas /(M gas + M„), where M gas includes Helium. The gas 
fraction represents the product of these opposing forces, and 



the evolution of / gas therefore encodes important informa- 
tion about their relative strength over time. 

Figure [9] shows the baryonic gas fraction for our sample 
of SMGs. To extend the redshift range of the plot we include 
four SMGs at high redshift [z > 4), as recently presented by 



Schinnerer et al. (20081, Daddi et al. (20091, Coppin et al 



(20101, and Riechers et al. (20101, in addition to GN20 in 



our sample at z — 4.05. We see a substantial variation in 
/ gas across the population, with values ranging from 0.1- 
0.9. However, overall we find the SMGs to be very gas rich, 
as expected, with a median f gas = 0.43T0.05. The evolution 
of the mean SMG / gas is overplotted on Fig. [9] We include 
on this plot the mean value for star-forming galaxies in the 
local Universe ( Bothwell et al.||2009 l as well measurements 
of UV-selected star- forming galaxies at z <; 1 from |Tacconi| 
et all (120101) andlDaddi et all d2010b, and mid-IR-selected 



z ~ 0.4 LIRGs from Geach et al. (2011). Note that these 
samples of galaxies use differing values of a appropriate for 
their respective types; Tacconi et al. (2010) adopt a — 3.2, 
Daddi et al. (2010) use a = 3.6, while Geach et al. (2011) 
and Bothwell et al. (2009) use a local "Galactic" value of 
4.6. 

It is notable that while SMGs are clearly highly gas- 
rich systems: at every epoch in our sample molecular gas 
represents 40-60% of their total baryonic content, they do 
not appear to be more gas rich than "normal" star-forming 
galaxies at comparable redshifts, despite being selected to be 
a strongly star-forming population. Indeed, although they 
have comparable baryonic gas fractions, the Tacconi et al. 
(2010) galaxies have mean SFRs of 100 M yr" 1 atz~ 1.2, 
and 150MQyr _1 at z ~ 2.3, lower than the SFRs of our 
SMG sample at the same redshifts (Table 5) by factors of ~ 4 
and ~ 8 respectively. The comparable gas fractions indicate 
that the SMGs are no less evolved - in terms of their gas 
properties - than less-active galaxies at similar epochs. 

It is also interesting to note the redshift evolution of 
the baryonic gas fraction. The median gas fraction of the 
sample remains approximately constant across the redshift 
range of our sample, down to z ~ 1.5. There is a steep de- 
crease in the mean gas fraction of galaxies below z ~ 1, 
however, with median gas fractions dropping to ~ 20% by 
z ~ 0.4 (Geach et al. 2011), and < 10% for typical star- 
forming galaxies in the local Universe (Bothwell et al. 2009). 
This drop below z ~ 1 coincides with both the observation- 
ally observed drop in global star-formation rate density (e.g. 
Hopkins fc Beacom|2006l |Sobral et al.|2011| - shown on Fig. 
9j), and the theoretically-predicted fall in specific accretion 
rates ( Krumholz fc Dekel|2011 1. It seems that with globally 
falling accretion rates below z ~ 2, the mean size of the gas 
reservoirs available to fuel starbursts decreases. 

In order to compare to simulations, we also include in 
Fig. [9] the predictions from the semi-analytic galaxy forma- 
tion model GALFORM ( Lagos et al.||2011 see also Swinbank 
et al. 2008). For this analysis, we applied a "LIRG" selec- 
tion criteria, so that only galaxies with Lfir > 10 Lq 
were considered. In order to have a reasonable comparison 
to the SMG population, we also only considered galaxies 
inhabitin g the most massive haloes, with masses of 10 12 - 
10 13 Mq ( Hickox et al|2012 l. It is clear that the real SMGs 
are somewhat gas deficient compared to those in the semi- 
analytical model, which predicts very high median gas frac- 
tions of > 70% above z ~ 2. Swinbank et al. (2008) drew a 
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similar conclusion from their comparison of GALFORM's pre- 
dictions for SMGs with observations: that the model had 
too many of the baryons in z ~ 2 sub-mm-selected ULIRGs 
in the form of gas, compared to their mass of stars. Thus 
it appears that the model needs to increase the efficiency 
of star-formation in massive galaxies at high redshifts to be 
match the observations. 

Other cosmological models fare somewhat better when 
compared to our observational data. The cosmological hy- 
drodynamic simulations of Dave et al (2011) predict smaller 
gas fractions for massive galaxies at high redshift - typi- 
cally ~ 30% for a galaxy with a stellar mass of 10 11 Mq at 
z ~ 2 closer to the current observations (the precise figure 
depends on input parameters, such as the choice of wind 
model). However, a major weakness of this comparison is 
that these model galaxies do not have a "high-SFR" cri- 
terion applied to them (to attempt to approximate SMGs) 
and so the relevance of this model comparison to the SMG 
population is not clear. 



4.6 A comparison of mass measurements 

As explained above, there are several properties of the SMG 
population which a detection survey such as ours is not ide- 
ally suited to study. These include the size of the 12 CO reser- 
voir, as well as the kinematic mode dominating the gas dy- 
namics (which we present as either rotational or virialised). 
These quantities can be studied in tandem, however, in order 
to try to explore the range of physically-motivated param- 
eters that describe the SMG population. As an example, 
it is unlikely that the SMG population as a whole is both 
rotation-dominated and compact]^] as the total dynamical 
masses derived under these assumptions are typically lower 
than our adopted stellar masses (which provide only one 
component of the total gravitational potential). 

It is possible to explore the range of parameters 
which result in physically-plausible outcomes. The dynami- 
cal masses (derived above) and the total baryonic masses 



require a value of a = 4.9 ± 0.8 - greater even than the 
value adopted for secular discs at z ~ 0, and inappropriate 



(our gas masses, plus stellar masses derived by Hainline 
|et al.|[20TT| of our SMGs effectively function as indepen- 
dent measurements of the total mass in the central regions 
of the galaxies (where the baryons are expected to domi- 
nate the kinematics; see j |4.1[ |. As a test of what constitutes 
a "physically plausible" set of parameters, we allow a to 
vary, calculating the value of a required in order to force 
agreement between the two independent measurements of 
the total mass. In other words, we assume 



1.36aL' co +M* 
M dy „ 



(12) 



(Where the factor 1.36 is to account for interstellar he- 
lium). Just as a combination of compact (3kpc) sizes and 
rotational kinematics results in a mean M* > Md yn (in our 
formalism above, implying a negative value of a), a com- 
bination of extended (7kpc) sizes and virialised kinemat- 
ics causes a similarly unphysical outcome. The dynamical 
masses derived in this circumstance are large enough to 



4 Which, as above, we take to mean a fiducial radius of 3 kpc; 
similarly, in this discussion we take "extended" to mean a fiducial 



for highly luminous starbursts, such as SMGs (see Solomon 
& Vanden Bout 2005 1 . Alternately, these large dynamical 



masses imply a dominant dark matter component to the 
kinematics of the SMG, something which is again unlikely 



radius of 7 kpc. 



(see p~ll. 

The two combinations that result in a physically moti- 
vated value of a, compatible (to within a factor of 2) to val- 
ues adopted typically adopted for ULIRGs/SMGs, are (1) 
compact sizes and virialised kinematics, and (2) extended 
sizes and rotational kinematics. Figure [To] shows the ratio 
of the total baryonic mass to dynamical mass for both of 
these models of our 12 CO-detected sample, plotted against 
redshift. 

There is a large scatter in each plot, indicative of the 
large uncertainties on each parameter (including the popu- 
lation averaged value of sin(i) adopted in the rotational es- 
timator). However, the scatter in the population is, in each 
case, larger than the expected la error (indicated in the bot- 
tom right), implying that there is a significant variation in 
the value of our derived mass ratio across the population. 
Given that the innate mass ratio is likely to be very close to 
unity (modulo a small dark matter contribution), this ad- 
ditional scatter, beyond the formal error, is likely a result 
of the globally-chosen parameters applying poorly to some 
individual SMGs - an incorrect value of a, for example, or 
an incorrect 12 CO radius. This suggests that the models be- 
ing discussed here almost certainly do not apply universally, 
across the entire SMG population - even in an average sense. 

As discussed above, it seems likely that the SMG pop- 
ulation is a heterogeneous mix of merger-like and disc-like 
galaxies, with the relative contribution of these to the total 
population varying as a function of luminosity. Considered 
within this framework, it could be argued that the SMG 
population contains both objects which are "compact and 
virialised" , and "extended and rotation-dominated" , with 
the former corresponding to the luminous "merger" sys- 
tems, and the latter representing the more "disc-like" ob- 
jects (such as the "Eyelash", SMM J2135— 0102). Solidifying 
this picture, however, will require higher-resolution measure- 
ments of 12 CO in a large number of SMGs, capable of re- 
solving the kinematic structure and the extent of the 12 CO 
reservoir. Using such measurements, it could also be possi- 
ble to use this technique to place dynamical constraints on 
the value of the 12 CO/H2 factor a for the SMG population. 



5 THE EVOLUTION OF SMBHS IN SMGS 

The characterisation of super-massive black holes (SMBHs) 
in high-redshift SMG hosts has the potential to shed light 
on the co-evolution of SMBHs and spheroids, an impor- 
tant facet of galaxy formation studies. It is well known that 
there exists a correlation between the mass of a supermas- 
sive black hole (Mbh), and the velocity dispersion (a) of 
its host spheroid: Mbh-o" ( Magorrian et al.|1998 1 , implying 
well-regulated feedback mechanisms in place which couple 
the two components. Indeed, the correlation is so good that 
the scatter is no more than would be expected from pure 
measurement error alone ( |Ferrarese fc Merr itt 2000). This 
tight correlation suggests a symbiosis between galaxy forma- 
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Figure 10. The ratio of total baryonic mass (the sum of the gas plus stellar masses) to dynamical mass, plotted against redshift, for 
both kinematical models of SMGs (rotational kinematics shown in the right panel, and virialised kinematics shown in the left panel). 
As the central regions of the SMGs — traced by the molecular gas - are expected to be baryon dominated we determine the value of the 
12 CO/H2 conversion factor (a) necessary to result of this ratio of unity (the values are reported in the respective panels). This figure 
also illustrate the scatter of the SMG population around this median ratio, which is larger than the l-cr error (shown in the bottom 
right), and may, in part, be driven by our assumption of a single model to describe to all SMGs. A median ratio of unity is indicated 
with a solid like, and dotted lines indicate a factor of 2 variation. 



tion and the formation of the central SMBHs, feedback from 
which is thought to play an integral part in the evolution of 
the most massive galaxies (e.g. Benson et al.|2003 l. 



5.1 AGN-dominated SMGs 



As reported by Alexander et al. ( 2005 1 , the high AGN frac- 



tion seen in the SMG population indicates relatively con- 
tinuous BH growth occurring throughout the intense star 
formation phase. This is in line with theoretical models of 



ULIRGs and SMGs (i.e. Narayanan et al. 2010 1: major merg- 
ers efficiently transport gas into the dense central regions, 
which - as well as fuelling a nuclear starburst - can effi- 
ciently "feed" a BH. 

While observational studies of SMGs have generally 
concluded that their prodigious bolometric luminosities are 



powered star-formation activity (e.g. Frayer et al.J 



Frayer et al. 



1998 



[19991 Alexander et al. 2005b; but see |Hill fc| 



Shanks 20111, many studies have found that a significant 
minority of SMGs do host a luminous AGN (Alexander et al. 
2005a, 2008 ; |Lutz et al.|2010| [Hamline et al.|2lml|Wardlowl 
et al.|2011| Ivison et al. 2011). There are a number of meth- 
ods that can be used to identify AGN activity, including 
X-ray emission ( Alexander et al.|2003 1 , mid-IR colour selec- 
tion ( Ivison et al. 2004 1 , mid-IR spectral properties ( |Valiante 
et al. 2007 Menendez-Delmestre et al. 20091, and optical 



spectral properties (Swinbank et al. 2004 1. Hainline et al 



(20111 analysed a sample of X-ray-observed SMGs, finding 



that the fractional contribution to the mid-IR from non- 
stellar (i.e. power-law) emission provides an excellent proxy 
for the hard X-ray luminosity and so provides a robust AGN 
diagnostic. 



We have compiled the AGN classifications for our sam- 
ple: we find 13/40 SMGs in our sample have been observed 
to have AGN activity, based on one (or more) of the above 
diagnostics. It must be noted, however, that the selection of 
the C05 sample could result in a bias towards SMGs con- 
taining AGN, as a result radio preselection, and/or the pres- 
ence of strong emission lines in their restframe UV spectra 
(this is discussed in |Hainline et aL]|2011| see also |Wardlow| 
et al. 2011 1. If this is the case, then our sample - being 
substantially drawn from the C05 sample could also be bi- 
ased towards a higher percentage of AGN than the general 
SMG population. Nevertheless, we can still use our sample 
to investigate the correlation between the presence of an 
AGN and its luminosity, with other physical properties of 
the SMG hosts. 

In fact, we find little or no correlation between the pres- 
ence of an actively-fueled SMBH, shown by AGN activity, 
and the 12 CO properties of the SMGs in our sample: gas 
masses, gas fractions, and FWHM-derived dynamical masses 
for the AGN-hosting-subsample are all consistent with be- 
ing identical to the 12 CO sample as a whole. This is in-line 
with the picture that SMGs are star-formation dominated, 
with activity from a central AGN playing little part in deter- 
mining the properties of the system. Indeed, it implies that 
the AGN activity occurs in all classes of SMGs, irrespective 
of dynamical or gas masses. The lack of any difference in 
the gas fractions of the SMGs showing AGN activity and 
those without, may be indicating that any feedback effect 
from the AGN influences both the dust and gas reservoirs 
in these galaxies. In this way those "evolved" systems where 
the AGN have had the most influence would no longer ap- 
pear in our 850-/xm-selected SMG sample due to an increase 
in dust temperature or a reduction in cold dust mass. 
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Figure 11. The Mbh~°"* relation for the SMGs in our sample 
with measured X-ray luminosities, which have been converted into 
SMBH masses following Alexander et al. (2005). We also plot 
the same quantities for local elliptical galaxies, as presented by 
Giiltekin et al. (2009), and the best fitting relation from the same 
work (with a slope /3 = 4.24). The SMGs have SMBH masses al- 
most an order of magnitude lower than their dynamical properties 
would suggest, suggesting that they represent a population yet to 
evolve on to the Mbh~ ct relation. The size of the SMG points 
scales with the baryonic gas fraction, two examples of which are 
shown in the key. There is no correlation between gas fraction 
and position on the diagram, suggesting that the rapid SMBH 
"feeding" stage may occur after the starburst episode. 



5.2 The Mbh-ct relation 

SMGs have been proposed as the precursors of the mas- 
sive ellipticals seen in the local Universe (e.g. Swinbank et 
al. 2006; Hickox et al. 2012). This population also exhibit 
the tightest Mbh-c relation (Giiltekin et al. 



2009) 



and as 

such examining their Mbh-c connection during the time of 
their peak star-formation activity may be useful if we are 
to understand the role of feedback in the galaxy formation 
process. Unfortunately, one difficulty with addressing the 
evolution of the Mbh-o" relation in SMGs using direct kine- 
matic measures is the fact that observations of the central 
stellar velocity dispersions at high redshift are not possible 
with the current generation of telescopes. Similarly, use of 
optical and near-IR emission line g tracer of host 

galaxy dynamics is complicated by the mixture of obscura- 
tion and potential outflow in these systems (e.g. Swinbank et 
al. 2006; Alexander et al. 2010). Previous attempts to study 
the Mbh-o" relation in SMGs have therefore been forced to 
use estimates of the stellar and gas masses as a proxy for a 
of the host galaxy (Borys et al. 2005; Alexander et al. 2005, 
2008). These studies uncovered an apparent offset between 
the relation for SMGs and that seen for spheroids at z ~ 0, 
in the sense that the SMBH masses were lower in the SMGs, 
relative to the host galaxy mass. This behaviour is the re- 
verse of that seen in studies of the SMBHs in QSO s jPeng 



et al.|2006[), including those using CO dynamics (Coppin 



et al. 12008), suggesting that SMGs must grow their SMBHs 



to lie on the present day relation. 
Our large sample of SMGs 



(providing well-constrained individual baryonic masses and 
robust kinematic information) is also a useful tool for 
analysing the relationship of the mass of the SMBH to 
that of their galaxy hosts, and compare this to the present- 
day Mbh-c relation to search for evolutionary changes. We 
therefore start by estimating the SMBH masses for the 
SMGs following [Alexander etH] p505) . We derive SMBH 
masses from the X-ray luminosity integrated from 0.5-8 keV 



(Alexander et al. 20031, assuming that the X-ray emission 
from the AGN accounts for 6^4 2 % of its total bolometric 
luminosity ( Elvis et al.||1994 l. Next, an "Eddington ratio", 
rj, can be adopted, which is the ratio of the AGN bolometric 
luminosity to the Eddington luminosity. As the Eddington 
luminosity is simply a function of black hole mass, this allows 
us to estimate the mass of the central SMBH. Of course, the 
adoption of an assumed Eddington ratio is a critical step 
here, and a source of uncertainty. Alexander et al. ( 2008 1 



investigated a sample of SMGs taken from the parent sam- 
ple of |Chapman et aL (20051, and estimate a typical value 
for the Eddington ratio of rj = 0.2-O.fj^] Here we adopt a 
value of rj = 0A J ^\ 1 and include the uncertainty in rj in the 
resulting uncertainty on the derived SMBH masses. 

We then derive the equivalent velocity dispersion of the 



galaxy from our 
Ho] ( [2Wfl ): 



CO linewidths, using the prescription of 



log 10 g = (1.26 ± 0.05) log {Jf-^j - (0.78 ± 0.11) 



(13) 



where a, is the stellar velocity dispersion of the bulge, 
and V20 is the velocity width at 20% of the peak. This 
can be somewhat uncertain for lower luminosity galaxies 
(and, of course, varies for sources with non-Gaussian emis- 
sion profiles), but for IR-bright galaxies (Lir ~ 10 12 Lq) 
such as our SMGs the relationship between V20 and FWHM 
(= V50) approaches the analytic relation for a Gaussian pro- 
file, V20 ^ 1-5 x FWHM ( PP007) . 



Figure 



11 



shows a versus Lx (and the resulting derived 
values of Mbh and Md yn ) for our sample of SMGs, which 
have been coded according to their baryonic gas fraction. In- 
cluded on the plot are local elliptical galaxies (thought to be 
the low-redshifts decedents of SMGs), as well as a fit to the 
local data (Giiltekin et al. 2009). The SMGs in our sample 
lie below the local Mbh-ct relationship, with SMBH masses 
approximately an order of magnitude lower than would be 
predicted based on their dynamical properties. This confirms 
the offset previously seen for SMGs using indirect estimates 



of their masses (e.g. Borys et al. 2005 Alexander et al. 2008). 

if we wished 



As discussed by Alexander et al 



with CO observations 



( |2008[ ), 

to remove this offset from the local Mbh-c by adopting a 
lower Eddington Ratio, we would be forced to an extremely 
low value of r\ ~ 0.05. This would imply an unfeasibly slow 
SMBH growth rate during a phase where the necessary gas 
supply should be abundant. We conclude that SMGs rep- 
resent a population which has yet to evolve onto the low- 
redshift Mbh-ct relation. While their high FWHMs (espe- 
cially in the high-J up lines, which trace the compact, cen- 
tral potentials) are indicative of massive galaxies, the cen- 



5 This range in estimations of r\ is primarily driven by the un- 
known distribution of gas within the Broad Line Region (BLR) 
of the black hole. 
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tral SMBHs are still relatively immature, having yet to grow 
substantially due to rapid mass accretion. 

An SMG's position in the Mbh-o" plane may evolve 
with age, as the central SMBH accretes mass and the gas 
is turned into stars. A low baryonic gas fraction could be 
interpreted as a sign that the SMG is a "more evolved" 
spheroid towards the end of the starburst episode, having 
already converted much of its gas reservoir into stars. As a 
result, if the SMBH growth occurs concomitantly with the 
starburst phase, we might expect gas fractions to correlate 
inversely with the ratio of SMBH to dynamical mass. As we 
see in Figure there is no apparent trend for lower / gas 
SMGs to be closer to the present-day Mbh-o" relation, al- 
though the small number of SMGs and the significant scat- 
ter which may be present in a due to geometrical projec- 
tion, weaken the strength of any conclusions. Equally, the 
lack of any apparent correlation may suggest that the rapid 
BH "feeding" stage occurs separately, after the starburst 
episode which causes the galaxies to be sub-mm selected. 
Indeed the timescales required to grow the SMBHs to the 
point that they resemble those o f their z ~ counte rparts 

al.|20 08[) are 
)right ] 



Alexander et 



(which have M B h > 4 x 10 s M q ; 

signficantly longer than the expected sub-mm-bright lifetime 
(100-300 Myr; |Swinbank et~aT1|2006] Hickox et al. 2012), 
again indicating the need for a subsequent phase of predom- 
inantly SMBH-growth, most likely associated with strong 
AGN activity. Interferometric observations of far-IR-bright 
QSOs could be a way to approach this phase from the "other 
side" , potentially shedding light on this interesting problem 
(Simpson et al. in prep). 



6 CONCLUSIONS 

We have presented results from a large PdBI survey for 
molecular gas in luminous SMGs. We observed 40 SMGs 
with well-defined redshifts, detecting 12 CO emission in 32. 
For the remaining eight SMGs without detected 12 CO emis- 
sion, we constrain the upper limits on their 12 CO flux. We 
used gas masses and dynamical masses from our sample, 
combined with ancillary multi-wavelength data (including 
far-IR luminosities and stellar masses), to discuss the physi- 
cal properties of the SMG population. Our main conclusions 
are as follows: 

• Analysing the median 12 CO spectral line energy 
distribution for our SMG sample, after normalisation to 
a mean far-IR luminosity, we find the SLED rises up to 

are sparse at J up > 5, but within the errors 
we find evidence for a turnover. 

• Plotting 12 CO luminosity against linewidth, we find 
that Lq O j 1 _ j oc a 2 is a good fit to the trend seen in our 
sample. The scatter around this relation is lower than would 
be expected for a population of randomly orientated discs. 
This suggests that the ISM in luminous SMGs typically 
exists as a thick disc or turbulent ellipsoid. 



respectively. In these calculations we have kept the depen- 
dence on the size of the 12 CO reservoir, which is not well 
constrained for the majority of our sources. We also find 
that 20-28% of the sample exhibit double-peaked 12 CO 
profiles, which we interpret as a signature of an on-going 
merger. 

• We use far-infrared luminosities to assess the star- 
formation efficiency in our SMGs, finding a steepening of 
the Lq -Lfir relation as a function of increasing 12 CO J up 
transition. While we find an approximately linear slope to 
the Lq O j 1 _ j-Lfir relation, without brightness temperature 
ratio measurements for individual sources this result is 
uncertain, and we suggest that it is difficult to draw any 
strong conclusions about the gas depletion timescales of the 
different populations from high-J up observations. 

• We see little evidence of evolution in the baryonic gas 
fraction in SMGs with redshift. The median value of our 
sample 40-60% at all redshifts. By comparison to recently 
derived gas fraction for other high-redshift populations, 
we conclude that SMGs do not have significantly higher 
gas fractions than more modestly star-forming galaxies at 
similar redshifts. 

• We compile a variety of archival multi-wavelength data 
in order to analyse the AGN contribution to the SMG 
population as a whole. We find that 30% of our sample 
have some indication of AGN activity, but the presence of 
an AGN does not seem to have a significant effect on the 
gas properties of the host SMG. We also make use of deep 
X-ray data to estimate SMBH masses for our SMGs. We 
find that the SMBHs in SMGs lie substantially below the 
z ~ Mbh-o" relation, and that there is no correlation 
between SMG gas fraction and SMBH mass. We conclude 
that the SMBH growth phase occurs separately from, and 
after, the "star formation" phase, which our sample of 
SMGs is undergoing. 
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• The 12 CO linewidths are, in general, broad, hav- 
ing a mean FWHM of 510 ± 80kms _1 . Using these 
line widths we derive dynamical masses and find mean 
values - for virial and rotational mass estimators - of 
(7.1 ± 1.0) x 10 10 iikpc M and (1.6 ± 0.3) x 10 10 R kpc M 
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Table 4. Observational parameters of the programme SMGs. (Full table available from authors upon re- 
quest). 

ID Transition CO position z I C o FWHM SgsOMmF] S£ 4GHz 

(J2000) (Jykms- 1 ) (km s" 1 ) (mjy) (/ijy) 

Detections 

SMMJ105141+571952 (2-1) 10 51 41.31 +57 19 52.0 1.2138 2.50 ±0.19 625 ± 50 4.6 295 



N on- Detection! 



SMM J105230+572209 (3-2) 2.6011 <0.23 11.0 86 



a From C05 (and references therein), unless otherwise stated. 

CO upper limits derived as detailed in §4.3 
c Redshifts for non-detections are quoted at the centre of the receiver bandwidth. 
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Table 5. Physical properties of the SMGs. (Full table available from authors upon request). 



ID 


T ' 

^CO 
(10 10 
K km s — 1 pc 2 ) 


L CO(l-oJj 
(10 10 
Kkms -1 pc 2 ) 


log 10 M(H 2 j^| 

(M ) 


l°gl0 M dyn 

(rotational) 
(^k P c M ) 


loglO M dyn 

(virial) 
(#kpc M ) 


log 10 M*g 
(M ) 


log 10 L PIR AGN^] 
(L©) 


Detections 

SMM J105141+571952 

N on- Detections 
SMM J105230+572209 


4.8 ±0.4 
<0.8 


5.8 ±0.5 
<1.5 


10.76 ±0.03 
< 10.17 


10.39 ±0.04 


11.04 ± 0.04 


11.34 


12.64 
12.90 



a Converted using an excitation model as detailed in §3 
b Assuming a = 1.0 M (Kkms -1 pc 2 ) -1 
c From Hainline et al. (2011) 

^ Letters represent AGN activity diagnosed using the following estimators: 
A = X-ray (Alexander et al. 2005) 

B = Ha (Swinbank et al. 2004; Menendez-Delmestre et al. 2009) 
C = Power-law NIR component fractions (Hainline et al. 2011) 
D = IRS spectral properties (Menendez-Delmestre et al. 2009) 
E = Colour selection (Ivison et al. 2004; Hainline et al. 2009) 
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Figure Al. SMM021725-0459. Note: All maps and spectra in the appendices are displayed identically to this, unless otherwise stated in 
the respective captions. Left panel: Multi-wavelength source map. Blue contours: integrated 12 CO image in la steps, starting at la. Red 
contours: 24/jm (where available). The background is a stacked image of all 4 Spitzer IRAC bands (where available; otherwise we have 
used other Near-IR imaging; A'-band, is shown here). The FWHM of the synthesised beam is shown in the corner, and the green cross 
indicates the centre of the radio emission. Right panel: 12 CO spectrum, taken at the point of maximum 12 CO flux, with the best fitting 
Gaussian profile overlaid in red. In the cases where a double-peaked profile provided a superior fit (as described in the text), the best 
fitting double Gaussian profile is shown. Vertical lines indicate previously obtained redshifts, derived as from the following wavelengths: 
black dashed line = from UV; red dot-dashed line = from Ho or other nebular emission line; green solid line = PAH. For the sources not 
detected in CO, the spectrum has been taken at phase centre. The SMM021725-0459 image has been stretched to account for the bright 
nearby star. 12 CO is detected here at 4.7a. The 850/xm source, SMM021725-0459 (Coppin et al. 2006), was followed up spectroscopically 
(Alaghband-Zadeh et al. 2012), detecting Ha as indicated, and the CO(4-3) and [CI] (1-0) spectra are published in (S. Alaghband-Zadeh 
et al. in prep). 
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Figure A2. SMM021738-0503. 12 CO is detected at 5.3 a. The 850^m source, SMM021738-0503 (Coppin et al. 2006), was followed up 
spcctroscopically (Alaghband-Zadeh et al. 2012) with Ha derived redshifts for both i^-band sources lying under the CO contours. The 
redshift of the northern A"-band source is indicated on the CO spectrum, the southern source lying +700 km/s offset. The CO(4-3) and 
[CI] (1-0) spectra are published in (S. Alaghband-Zadeh et al. in prep). 




Figure A3. SMM021738-0505. 12 CO is detected at 6.2 a. The 850Atm source, SMM021738-0505 (Coppin et al. 2006), was followed up 
spectroscopically (Alaghband-Zadeh et al. 2012), detecting Ha as indicated, and the CO(4-3) and [CI](l-0) spectra are published in (S. 
Alaghband-Zadeh et al. in prep). 
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Figure A4. SMM030227+0006. 12 CO is detected at 4.2 a, and is lined up well with both IRAC and 24/xm emission. The PAH redshift 
is slightly inconsistent, but consistent with dz ~ 0.02 redshift fitting errors (Menendez-Delmestre et al. 2007, 2009) (z = 1.408 from UV, 
z = 1.4076 from Ha, z = 1.43 from PAH). 
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Figure A5. SMM044315+0210. 12 CO is detected at 7.0 a. This SMG was previously published by Neri et al. (2003); Greve et al. (2005); 
Tacconi et al. (2006). 
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Figure A6. SMM094304+4700. 12 CO is detected at 11.0 a. This SMG was previously published by Neri et al. (2003); Greve et al. 
(2005); Tacconi et al. (2006). In Engel et al. (2010), the H6 component of this 30 kpc separated system is also detected at 5cr in CO (6-5). 




Figure A7. SMM105141+5719. 12 CO is detected at 10.7cr, and CO contours are shown double spaced (2<r, 4a, 6<r, etc). The 24/^m 
emission is well aligned with the CO. In addition, high spatial resolution (A-config) CO(4-3) data for this SMG is published in Engel et 
al. (2010), resolving the CO into two components. 
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Figure A8. SMM105151+5726. 12 CO is detected at 6.2a, but near the edge of the bandpass. All emission components are aligned to 
the radio centre within errors. The published UV-derived redshift, z = 1.147 in Ivison et al. (2005) results from a foreground UV-bright 
galaxy. The PAH redshift can be difficult to constrain in this redshift range, and we show the best two fits - the higher of the two 
estimates was the initial estimate from Menendez-Delmestre et al. (2009), while the lower estimate is in agreement with the CO redshift. 
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Figure A9. SMM105227+5725. CO is detected at 5. Oct, though it is ~5" north of other emission components. 
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Figure A10. SMM105307+5724. 12 CO is well detected at 6.9a, and all emission components are close to radio centre. In addition, high 
spatial resolution (A-config) CO(3-2) data for this SMG is published in Bothwell et al. (2010), resolving the CO into a north-south, 1.5" 
extended source. 




Figure All. SMM123549+6215. CO contours are double spaced (2a, 4a, 6a, etc). 12 CO is well detected at 9.3a. The PAH redshift is 
inconsistent with other determinations (2 = 2.203 from UV, 2 = 2.2032 from Ho; compare to 2 = 2.24 from PAH), slightly larger than 
the dz = 0.02 rms error on the PAH fit. This SMG was previously published by Tacconi et al. (2006). 
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Figure A12. SMM123555+6209. 12 CO is detected at 4.3cr. The peak of emission is between 6 to 8 arcsec to the south of the 
IRAC/radio/24/xm positions. Both the published Lya-based UV redshift and the PAH-derived redshift are higher than the CO: z = 1.875 
from UV and z = 1.88 from PAH, while the 12 CO line is detected at z = 1.864. Followup spectroscopic analysis suggests a slightly lower 
systemic redshift z = 1.868 based on weak inter-stellar absorption features blueshifted relative to the Lyo emission line. 
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Figure A13. SMM123606+6210. 12 CO is reasonably well detected at 4.3(T (approximately 5" to the NE of the radio centre) but there 
is no apparent near-IR counterpart. A radio source is detected by Spitzer-MIPS and IRAC, 20 arcseconds away to the NE, but is likely 
unconnected with the SMG. 
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Figure A14. SMM123618+6210. 12 CO is detected at 4. Oct close to the radio centre and within the elongated beam resulting from the 
poor UV-plane coverage, here marked marked with a purple star. The green cross indicates here the 850/^m centre (which was chosen 
for pointing). We classify this as a candidate detection as in the text, as a result of the offset, CO significance, and lack of an optical 
redshift to verify the PAH. 
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Figure A15. SMM123618+6215. CO contours are double spaced (2a, 4cr, da, etc). 12 CO is strongly detected at 8.6rr and aligned with 
the 24/xm emission. As discussed in Chapman et al. (2009) and Bothwell et al. (2010 - where higher resolution CO(4-3) data is published), 
the redshift originally derived from UV is for the galaxy 0.8" to the west, which led to a CO non-detection in a separate observation. 
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Figure A16. SMM123632+6207. 12 CO emission is detected at 4.5cr, offset by ~ 7" from the radio centre, but still consistent at 2a 
within the elongated CO beam. 




Figure A17. SMM123634+6212. CO contours are double spaced (2a, 4a, 6a, etc). 12 CO is strongly de tected at 9.5a. All emission well 
aligned with the radio counterpart. This source was also observed in 12 CO (2-1) by |Frayer et al.| ( |2008^ , who report a 12 CO luminosity 
and FWHM that matches the values reported here. Higher resolution CO (6-5) data was also published in Engel et al. (2010), where they 
marginally resolve a velocity field within the essentially unresolved source (HWHP 2.2 kpc). 
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Figure A18. SMM123707+6214. 12 CO is well detected at 5.9a. All emission components line up well. High resolution imaging (Tacconi 
et al. 2006, 2008) show this source to consist of two 12 CO emission components, separated by ~ 2.5", mirrored by the double radio and 
IRAC sources. 
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Figure A19. SMM123711+6213. The CO(3-2) detection for this source was originally published by |Chapman et al.|p008^ . The spectrum 
shown is the NW 12 CO source. 1 CO is detected at ~ 4.8cr - well aligned with radio, 24/xm and IRAC emission. The green cross here 
indicates the phase centre lying between the two 12 CO sources (see SMM123712+6213 below). SMM123711+6213 is also well detected 
at higher resolution (A-config) in 12 CO(4 — 3) at > 6cr, as reported in Bothwell et al. (2010), with lower resolution D-config 12 CO(4 — 3) 
data presented in Casey et al. (2011). 
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Figure A20. SMM123711+6222 = GN20. The phase centre (green cross) was chosen to lie between GN20 (east), and GN20.2 (west; 
not detected). Both galaxies are detected in lower-J imaging (Daddi et al. 2009). The spectrum consists of two lines; 12 CO(7-6) (fitted 
with a Gaussian function), and Cl, approximately 700 km/s bluewards of this. The data was previously presented in Casey et al. (2009), 
although their overestimated 2mm continuum strength led them to question the detection of the CO(7-6) line. 




Figure A21. SMM123712+6213. The spectrum shown is the SE 12 CO source. 12 CO is detected at ~ 4.7ct, but the source is offset from 
5" to the south of the radio position, and the mid-IR IRAC position. However, this source is also detected at 4.5cr in 12 CO(4 — 3) at the 
same southern offset position, as reported in Bothwell et al. (2010), lending additional weight to the reality of the 5" offset. 
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Figure A22. SMM131201+4242. 12 CO is detected at ~ 5.7<r. The image is taken from Greve et al. (2005). Higher resolution CO(6-5) 
data is published in Engel et al. (2010), resolving the CO into three components. 
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Figure A23. SMM131208+4241. The source is weakly detected in 12 CO at 3.6<r, but is aligned with the radio counterpart, and with 
the weak IRAC identification where the optical redshifts were measured. We identify this as a candidate detection mainly based on the 
weak CO detection. 
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Figure A24. SMM131232+4239. 12 CO is detected at 4.9<r, but is offset to the north by 8". The Ha redshift (obtained for the radio- 
identified IRAC source shown) is significantly offset from the CO as well — (z = 2.300 from Ha, while z = 2.320 is measured from UV 
Lya line, and the 12 CO line gives z = 2.332.). This source could represent a spatially extended merging system, similar to the SMG 
complex reported by Ivison et al. (2011). We identify this as a candidate detection based on the large offsets spatially and in velocity. 




Figure A25. SMM163639+4056. 12 CO is detected at 4.5cr, is offset 6" to the east of the radio position, but is consistent within 2a 
within the elongated beam shape. We identify this as a candidate detection. 
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Figure A26. SMM163650+4057. 12 CO is strongly detected at 8.1ct. The source is lined up with 24/xm and Near-IR emission at the 
position of the radio counterpart. This is a higher resolution image, and as a result the map size is shown as 25x25". The data was 
previously published by Tacconi et al. (2006, 2008) along with higher resolution CO(7-6) data showing an elongated linear feature. 




Figure A27. SMM163658+4057. 12 CO is detected at 6.3<x. The PAH redshift is slightly high, at z = 1.20 (but still consistent, given 
the typical dz ~ 0.02 error in fitting PAH features). 
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Figure A28. SMM163658+4105. 12 CO is well detected at 5.0<t. Source is well aligned with 24/xm emission. This source was previously 
published by Greve et al. (2005) and Tacconi et al. (2006, 2008), the higher resolution, higher J data showing a very compact CO(7-6) 
source (HWHP=0.8 kpc). 
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Figure A29. SMM163706+4053. The background image is / band, rather than a composite IRAC image. 12 CO is detected at 4.4er and 
aligned with the radio counterpart. The source was previously published in Greve et al. (2005). 
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Figure A30. SMM221704+0021. 12 CO detected at 6.0a, and is lined up with the emission, and two IRAC sources, both of which 

are detected in Ho. The PAH redshift (z=2.55) is slightly inconsistent (Ha and UV give a similar redshift to the CO, z = 2.517), but 
again within the typical PAH fitting errors. This source is detected in 12 CO (4—3) and in [CI](l-0) in Alaghband-Zadeh et al. (in prep). 




Figure A31. SMM221735+0015. 12 CO is detected at 4.7<x. There is no 24 (im emission associated with this source (Menendez-Delmestre 
et al. 2009). The PAH redshift (z = 3.21) is inconsistent with the CO data, however the UV spectrum implies a redshift consistent with 
the CO, z = 3.098. The source was previously published in Greve et al. (2005). 
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Figure A32. SMM221737+0010. The background image is if -band. 12 CO is detected at 4.2ct, approximately 5" to the east of the radio 
centre. There is extended 24um over the CO region, and a A-band source 3" to the West of the CO emission, from which the optical 
and near-IR spectra were obtained. We identify this as a candidate detection based on the weak CO detection and the offset from the 
optical source. 
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Figure Bl. SMM105230+5722. Non-detection. All the sources shown in Appendix B are true non-detections, which allow useful con- 
straints to be placed on the CO flux. 




Figure B2. SMM123600+6210. Non-detection. 
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Figure B3. SMM123621+6217. Non-detection. 



12 36 29.171 +62 10 45.5oj; 



o 



o 



30 20 10 -10 -20 -30 

RA (arcsec) 



1.0063 

6^ 



■At 
-1000 



Redshift C0( 2 - 1 ) 
1.0108 1.0153 



1.0197 



-500 

Velocity (km/s) 



500 



1000 



Figure B4. SMM123629+6210. Non-detection. 
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Figure B5. SMM123712+ 6212. Non-detection. 
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Figure B6. SMM163631+4055. Non-detection. 
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Figure B7. SMM163655+4059. Non-detection. The emission to the SE of phase centre is most likely spurious, as indicated by IRAC 
photometric redshift estimates. 
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